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It has been shown recently that intermittency of the Gledzer Ohkitani Yamada (GOY) shell 
model of turbulence has to be related to singular structures whose dynamics in the inertial range 
includes interactions with a background of fluctuations. In this paper we propose a statistical theory 
of these objects by modelling the incoherent background as a Gaussian white-noise forcing of small 
strength T. A general scheme is developed for constructing instantons in spatially discrete dynamical 
systems and the Cramer function governing the probability distribution of effective singularities of 
exponent z is computed up to first order in a semiclassical expansion in powers of F. The resulting 
predictions are compared with the statistics of coherent structures deduced from full simulations of 
the GOY model at very high Reynolds numbers. 
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I. INTRODUCTION 

(N ■ 

Are structures (sheets or filaments of vorticity) a vital ingredient of intermittency in 3D-incompressible turbulence? 
To date, this important question remains open M, and an answer starting from first principle, i.e., from a controlled 
approximation to the Navier-Stokes equations, seems over the horizon. The new understanding of the anomalous 
scaling in the Kraichnan's model of passive advection || , based on the identification of zero modes in the homogeneous 
Hopf equations for equal-time correlators, has in particular strengthened the belief that field-theoretical methods would 
Q\ j eventually be able to capture the full statistics of turbulent flows without an explicit account of structures. 

■ Interestingly enough, the relative interplay between coherent ordered structures and incoherent turbulent fluctua- 
tions turns out to be a subtle matter already in the restricted framework of the so-called shell models of turbulence 
||. It was noticed very soon [Q that elementary bricks of intermittency in those deterministic ID-cascade models 
could be pulses or bursts of activity growing in an almost self-similar way as they move from large to small scales. 
However, genuine dynamically stable self-similar solutions of the equations of motion in the inertial range display 
an unique scaling exponent (to be denoted below as zo), provided they are localized in k-space (which, in the shell 
model approach, reduces to a discrete set of wave numbers k n — 2^ n ~ x \ where the shell index n goes from 1 to oo). 
Furthermore, the exponent zq, giving the logarithmic slope of the velocity gradient spectrum left in the trail of the 
pulse, happens to be rather close to the Kolmogorov value 2/3 (zq — 0.72) in the case of the Gledzer-Okhitani- Yamada 
(GOY) model, in the range of parameters where it reproduces at best the multiscaling properties of real turbulent 
' flows. 

In Ref. H , the role played by the interaction of pulses with the rest of the flow in producing more singular events 
was unravelled, and a two-fluid picture was introduced, where coherent structures form in and propagate into a 
featureless random background. Our goal in this paper is to elevate this still rather qualitative proposal to the rank 
of a semi-quantitative theory and to test its predictive power about intermittency in the GOY model. We shall 
assume that turbulent fluctuations on the shells downstream the pulse, i.e., small scales, act on the coherent part of 
the flow as a random, white-in-time, Gaussian forcing and ask whether the inviscid stochastic extension of the GOY 
model obtained in this way is able to reproduce the statistics of strong deviations of the full turbulent system in the 
inertial range. There is a priori quite a lot of freedom in the parameterisation of the forcing. Therefore, in order to 
keep things as simple as possible, we bind ourselves to use a single adjustable parameter (hereafter noted T), which 
measures the level of noise. We consider the semiclassical limit T <C 1 of these systems and study the statistics of 
singular structures appearing in this regime. 

Semiclassical (or instanton) techniques are well suited to capture large and rare excursions of fluctuating fields . 
As such, they have gained recently a renewal of interest in the field of turbulence and have already led to noteworthy 
results in the context of Burger's turbulence and of the Kraichnan's model of passive scalar advection |ll|,^2[ . 

One usually starts from a path integral representation of high order structure functions and uses a saddle point 
approximation to determine the coupled field-force configurations contributing mostly to those quantities. The nature 
of the statistical object to be computed imposes precise boundary conditions on the physical field and the random 
force (respectively at large and small scales, where the cascade processes start and end). Instantons, which in the 
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inertial range often reduce to a self-similar collapse along some spatial dimensions, are eventually selected by a delicate 
matching procedure at the two boundaries. In shell models we are dealing with an intrinsically discrete lattice of 
logarithmic scales. As a consequence, the analytic computation of instantons is completely out of reach in the inertial 
range, not to speak about the matching on both sides of the cascade. To circumvent this difficulty, we shall focus 
on the probability distribution function (pdf) of scaling exponents after n cascade steps, P n (z), and argue that, in 
the semiclassical limit V <C 1, this pdf builds up from the neighborhood of a single self-similar instanton (of scaling 
exponent z) that dynamic stability considerations will help us to construct numerically. In order to get non trivial 
physics, it turns out to be necessary to perform the semiclassical expansion of — linin^oo — In P n (z) (the rate of 
rarefaction of singularities of scaling exponent z in the multifractal picture) up to next to leading order in powers 
of r. This can be achieved via a summation over quadratic fluctuations around the instantons, once the proper set 
of boundary conditions for the corresponding trajectories in configuration space has been defined. We shall show in 
details how to carry out this program and end up with a prediction for P n (z) lending itself to a straight confrontation 
with the pdf of effective scaling exponents of coherents events that can be extracted from simulations of the GOY 
model at very high Reynolds numbers. Although our interest lies primarily in gaining a better understanding of 
intermittency in the framework of shell models of turbulence, the emphasis will be put in this paper on the technical 
aspects of the method that we had to develop for computing instantons. We believe that this method is general 
enough to find applications in other contexts or physical problems, like for instance the motion of complex objects or 
excitations on ID-lattices in the presence of a co-moving random environment. 

The paper is organized as follows. In Section 0, we define the stochastic extensions of the GOY model that we 
shall study. In Section III the equations of motion for instantons will be derived using the well-known Martin Siggia 
Rose path integral representation of probability distribution functions for stochastic dynamical systems. Section IV 
is devoted to the computation o f self- similar extremal trajectories, with the theoretical considerations underlying the 
solution explained in Subsection IV A and its practical implementing, together with the results, exposed in Subsection 



IV B. The important effect of quadratic fluctuations and the rather heavy formal work behind their computation are 
discussed in Section [y| The comparison of the results issuing from the instanton approach with numerical data on 
the statistics of coherent structures in the genuine GOY model is given in Section [Vj. We conclude in Section VII , in 
particular as to the relevance of a two-fluid description of intermittency in shell models of turbulence. 



II. DEFINITION OF THE STOCHASTIC DYNAMICAL SYSTEM 



Equations of motion for the GOY model in the inertial range read : 

Q 2 (i - z)K-2<-i + <-iK+i - Q~ 2 K+iK+2, 



db n 
~~dt 



(2.1) 



where the complex variable b n — k n u n should be understood as the Fourier component of the gradient velocity field 
at wavenumber k n — Q n and the integer n runs from to +oo. Throughout this paper, usual values of parameters 
e = 0.5 and Q — 2 will be assumed. It is convenient to cast (2.1) in a vectorial form 



db 
~dl 



N[b], 



(2.2) 



where the infinite-dimensional vect or b is built up from the 6 n 's, while the n component of the nonlinear kernel N[b] 
is given by the right hand side of (2.1). It is worth noting at this point that b* ■ b = J2^=o l^n| 2 plays dimensionally 
the role of enstrophy in real flow and that the inverse square root of this quantity sets the order of magnitude of the 
smallest time scale on the shell lattice. 

Since quadratic nonlinearities lead generically to finite time singularities, it is very useful to introduce a desingu- 
larizing time variable r related to the physical time t by the differential law 



This turns (2.2) into 



db _ N[b] 

~ (b* -b) 1 ^' 



(2.3) 



(2.4) 
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where both sides of the equation have the same scaling dimension in the field b, which shows that an infinite "time" 
is now required to form a singularity by travelling across the whole shell axis. 

^From previous work Jl3[ ] , we know that every initial condition of finite enstrophy, when evolving under dynamics 



(2.4), eventually organizes itself in a soliton-like pulse, moving from large to small scales at a constant speed with 



an exponential growth of its amplitude. The asymptotic state is unique, up to trivial phase symmetries of the GOY 



model jlj], time translations or multiplicative rescaling of the field b, which all leave the equation of motion (|2.4| ) 
invariant. We may restrict our attention without loss of generality to the case where the phase pattern along the shell 
axis does not break into a three-sublattice structure. The asymptotic Floquet state, to be noted henceforth b°(r), is 
then purely real and such that : 

b n+1 (r + To) = exp(A T )b° n (T). (2.5) 

The period To is the "time" needed for the center of the pulse to go from shell n to shell n + 1, while the (positive) 
Lyapunov exponent Ao controls its growth. Both quantities To and Aq, are dynamically selected in an unique way. 
The scaling exponent zo associated to the pulse (fixing in particular the logarithmic slope of the spectrum left in its 
trail) can be extracted from the identity Q z ° — exp(A T ). Its value turns out to be 0.72 in the case of the GOY 
model for the choice of parameters stated before. 

We turn now to the stochastic models, that we are interested in solving by the instanton method. Their physical 



we assume that pulses para mete rize adequately singular (and temporally 



motivation has been explained in Ref . 
coherent) structures in shell models but that the deterministic dynamics (2.2) should be enlarged towards a stochastic 
one, in order to describe the interaction of a given pulse with incoherent fluctuations produced by the relaxation of 
the trails left by its predecessors. We are therefore led to consider the following extension of the original inviscid GOY 
model : 

^ =N[b] +VT(b* •b) 3 / 4 B[C]7 ? , (2.6) 
at 

where r\ is a Gaussian noise, delta-correlated in time and shell index, whose correlations read : 

«(*)»*.' (0) =Snn'S(t-t'). (2.7) 



The various factors coming in front of rj in ( p.6[ ) have the following meaning : the number T fixes the relative strength 

of incoherent fluctuations with respect to coherent ones and we shall be interested in the semiclassical limit of small 

r amenable to semi-analytic treatment. As will be clearer in a while, the overall scale factor (b* • b) 3 / 4 is there to 

keep noise relevant all along the cascade, thereby preserving scale invariance. Finally the matrix B[C], of zero scaling 

dimension in the field b since it depends only on the unit vector C = , b , may be used either to introduce spatial 

V( b *' b ) 

correlations of noise (along the shell axis) or to localize its action with respect to the instantaneous position of the 
pulse. Although the formalism to be developed in this paper can deal with the most general situation, we restricted 
ourselves in practical investigations to diagonal matrices S[C], just playing with the degree of localization of noise. 
Results will be presented for three rather emblematic choices of B : (i) B nn = 1, which describes a completely 
delocalized noise; (ii) B nn = C , *_ 2 C I *-i which keeps some flavour of the original GOY dynamics and makes noise 
active just at the leading edge of the pulse; and finally (iii) B nn = | C rl _ 5 1 2 + |C„_4| 2 , which removes the action of 
noise further away from the center of the pulse. We must emphasize that these particular choices were not dictated 
by rigorous considerations on the underlying dynamics of the GOY model, but rather used to scan the variety of 
behaviours which may be expected from such stochastic dynamical systems. It should be noted that the strucutre of 
the matrix B is not constrained by any conservation law, since the coherent part of the flow does not form anymore 
a closed system, even in the inertial range, in our two-fl uid d escription. Finally, to simplify the following analysis, 



we are going to restrict the fields b and r] in (2J3) and (2/7) to being real- valued vectors and neglect the effect of 



imaginary fluctuations. This is certainly not a serious restriction as for the instantons themselves, which are expected 
to be, like the self-similar deterministic solution described above, purely real, up to trivial phase symmetries of the 
GOY model. It can also be remarked that the model (ii) (which will be found later on to give the more convincing 
results) does not require a complex noise, since the phase degrees of freedom have already been incorporated in the 
definition of the matrix B in that case. 

While the deterministic dynamics (|2.2|) s elects a single self-similar solution blowing up in finite time with scaling 
exponent zq, the presence of noise in ( |2.6| ) allows for a continuum of scaling exponents, even in the manifold of 
normalizable fields b. In the small noise (or semiclassical) limit r « 1, the probability density of developing an 
effective growth exponent z after n 3> 1 cascade steps will take the form : 
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P n {z)~y/n exp 



sq(z) 

r 



si(z) 



(2.8) 



where Sq(z) is the action per unit cascade step of the self-similar extremalsolution of scaling exponent z of optimal 
bare Gaussian weight (or instanton), and si(.z) measures, to lowest order in T, how the basin of attraction of the 
instanton in phase space evolves with the number of cascade steps. Note that the argument of the exponential in 
fl2.8p , —(^Y^- + si(z)), is nothing but the Cramer function introduced in the theory of large deviations, which governs 
the rate of rarefaction of singularities in the multifractal picture [p"5f . We will show in this paper how to compute 
in a clean way both quantities sq(z) and s\(z) . B efore doing this, we must carefully handle problems related to the 



time discretization of the stochastic equation (2.6) since a consistent treatment of them is necessary to get the right 
expression of the first correction si(z). We shall adopt the view that the initial stochastic equation (2.6) is to be 
understood in the Stratonovich sense |16| . However, in the path integral formulation of stochastic dynamical systems 
that we shall heavily use in the followi ng, it is much simpler to work with the Ito prescription which, in the limit of 
small time steps, amounts integrating (^_^) within a basic Euler scheme with all b-dependent quantities in the r.h.s. 
estimated at the prepoint. When switching to the Ito discretization recipe, the stochastic equation has to be changed 
into : 



Vf(b-b) 3 / 4 B[C]r/, 



(2.9) 



where the new kernel Nr[b] differs from N[b] by the addition of the so-called Ito drift-term. We give, for the sake of 
completeness, the resulting expression of the n th component of Nr : 



N rn [h] =N n [h] 



lrA [(b .b) 3 /%](b.b) 3 /% 



(2.10) 



At this point, we may write down the discrete analog of ( |2.3| ) as At; = r i+1 — Tj = (bj • bj) ' (ij+i — U) (where i is the 
ti me in dex) and redefine the noise as 77 ^ — > = VT(hi ■ b^) -1 / 4 ^. This leads to the following stochastic extension 



(2.11) 



of (|2.4|) which will be the starting point of our formal analysis 



dh 



N r [b] 

(b-b)!/2 



{h-h) l i 2 B[C]d. 



with 



(2.12) 



III. EXTREMAL TRAJECTORIES FROM PATH INTEGRAL FORMULATION 

Statistics of classical fields in the presence of random forces can be examined with the help of field theoretical 

techniques formulated in fl7|| . In particular, the probability to go from point bj n at time r = to point by at time 
Tf may be written as a path integral : 

P(b in ,0;bf,T f ) = JvbVp exp-5[b,p], (3.1) 

where S[h, p] is an effective action to be defined below, p an auxiliary field conjugated to the physical one b and 
VhVp stands for : 

11 T^vT- [6 - Z) 



(2tt)<* 1 = 1 (2tt)<* ' 

In the last equation, the time interval Tf was divided into N subintervals of length At = 2£ (with h in = b and 
hf = bjv) and the num ber of shells was set to a finite value d, in order to give a clear meaning to the measure. For 



the problem of interest (2.11), the effective action S takes the form : 

. N r [b t ] \ : 
(b.-b,) 1 / 2 ; 2 



S[b,p]= »P<- ( b '+ 1 b ' - At ( u. .h.u/2 ) + 9 ( b « ■ b ') P' ■ B ^ tB M^ ■ ( 3 - 3 ) 



i=0 
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or, in the continuum limit 



r T f 

5[b, p] = / drip. 
Jo 



dh 

d^ 



N r [b] 
(b-b) 1 ^ 



: (bb)pS[b]*S[b]p. 



(3.4) 



The last term in ( |3.4| ), quadratic in p, appears as a result of averaging over the Gaussian noise £, while the first one, 
linear in p, would still be there in the absence of noise as a formal way of enforcing the deterministic equation of 
motion of b. 5[b, p] will be referred to in the following as the Martin-Siggia-Rose (MSR) action. 

Rescaling the auxiliary field p as p'/r puts an overall large fac tor 1/T in front of the effective action and opens 
the way to a saddle point approximation to the path integral (3.1). Extremization of the action with respect to the 
configurations of both fields b and p between times and r/, for fixed endpoints, leads in a straightforward way to 
the following set of coupled equations defining extremal trajectories : 



(b • b)B l B6 . 



dh _ N r [b] 
d^ ~ (b-b)V2 

dO 1 

^ = -*M - -0b[(b • b) 6 ■ B *B6\ 

UT 2 



In the above equation, we set p' 
<9 b N r [b] N r [b](g>b 



-iO and M. is the Jacobian matrix of the kernel 



N r [b] 
(b-b)V 2 



(b • b)V2 (b • b) 3 / 2 



As usual, equations (3J5) and (|3.6| ) inherit a canonical structure 



dh _ dH 
d^ ~ ~d0 



where the Hamiltonian TL reads 



H = 



dO 

dr 

0-N r [b] , 1 

(b ■ b)!/2 ' 



on 

dh ' 



^(b-b)(e -B*Be) 



(3.5) 
(3.6) 
M = 

(3.7) 
(3.8) 

(3.9) 



Since H. is not explicitly time-dependent, we conclude that its value, to be called the pseudo-energy in the sequel, is 
conserved along any extremal trajectory. The action S\b, 9] may be rewritten in terms of H. as 



s[b,e] 



dr e 



dh 
'~dr 



H 



(3.10) 



from which it is seen that the further requirement that the trajectory be extremal with respect to time reparametriza- 
tion leads to the condition of vanishing pseudo-energy 7i = 0. Noting that each term of TL in (3.9) has the same 
scaling dimension in b and (cither 1 or 2), one finds that 



dr 



(3.H) 



which shows that the overlap h-9 between the physical and auxiliary fields is also conserved, together with Tt, along an 
extremal trajectory. This property reflects the scaling invariance of the stochastic cascade processes we have in mind. 
We should at this point insist on the fact that, in contrast to instantons in the framework of equilibrium statistical 
mechanics or quantum mechanics, equations for extremal trajectories in stochastic dynamical systems describ e the 
real motion of the physical field in a particular "optimal" realization of the noise. The comparison of equations (2.11) 
and (3.5) shows indeed that the following relation holds between £ and 9 : 



£= (b-b) 1/2 *B6>. 



Like their deterministic parent ( |2.4| ), the equations of motion (3.5) and (j; 
solutions, such that 

b n +i(r + T) = expAT6„(r) , 
9 n+l {T + T) = exp-AT6 n ( T ), 



(3.12) 

j) sustain formally traveling wave-like 

(3.13) 
(3.14) 



5 



whose scaling exponent z = is expected to be now related to the overlap fi\ —b-8 (with z = zq for fii — 0, in the 
absence of noise). However there is little hope to fin d th ese solutions by a direct forward in time integration of (3.5) 
and (|3~6|), as could be done successfully for equation (2.4). This is b ecause the auxiliary field 6 intrinsically propagates 
"backward" in time, as is clear from the discretized version of ( |3.6| ) (deduced from the extremization of ( |3.3| )). I n the 
present problem, we have observed numerically that regular Floquet states emerge as dynamical attractors of (3.5) 
and ( |3.6| ) only for rather high values of /ii (otherwise the system evolves in a chaotic manner). They form a branch 
of solutions definitely distinct from the one to be obtained in the next Section and correspond presumably to local 
maxima of the action rather than the local minima of interest to us. 



IV. AN ITERATIVE METHOD FOR COMPUTING SELF-SIMILAR INSTANTONS 



A. Theory 

The previous considerations suggest that equations (3J)) and (3J:) should not be treated on the same footing. 
The careful examination of physical properties that instantons sho uld p ossess will give us keys for computing them. 
Assume for a while that a solution has been found, obeying to (|3. 13|) and (|3.14[) . We note b°(r) and £°(r) the 
corresponding configurations of b and The linearization of the equation of motion ( |3.5| ) at fixed noise leads to the 
following evolution of fluctuations 6h of b around b° : 



-^-<*b = C8b=M8b + (Sh ■ d b ) ( (b • b) 1 / 2 B£° 
dr V 



b" 



(4.1) 



The periodicity properties of the linear operator C e nsur e that the fluctuations of b may be decomposed on a complete 
set of eigendirections *j r (r) evolving according to (4.1) and such that 



*i r (r + T) 



<7i T 



T+i *„(r), 



(4.2) 



where T + i denotes translation by one unit in the right direction along the shell lattice. In practice we shall have to 
work with a finite number of shells d and, in order to get rid of boundary effects, it will be necessary to fully periodize 
the shell lattice : the index i then runs between 1 and d and the translation operator is easy to represent as a matrix. 
Formally, the SI^'s can be determined at time r = by diagonalizing the Floquet operator : 



Ut = 71 1 cxp / £dr, 



(4.3) 



where exp is a chronologically time ordered product (initial time on the right). One observes that b° satisfies ( |4.2| ) 
with a time averaged Lyapunov exponent a — A. 

We claim now that every initial condition b(0) evolving in the configuration of noise £ should be attracted towards 
the instantonic trajectory. If it were not true, some perturbations would be able to grow in the comoving frame of the 
pulse, thereby generating scaling exponents larger than z at no cost in the action, in contradiction with the hypothesis 
that the optimal realization of a singularity of exponent z has been found. This strong criterion of dynamic stability 
is another way of stating that the Cramer function should be insensitive to the details of the production of pulses in 
the forcing range. It implies that A is an upper bound for the real part of the oVs. Arranging the eigendirections 
in order of decreasing Re<7i, we are therefore led to identify b°(r) with ^i r (r). In the case of zero noise where we 
recover the deterministic solution of Section [n] (with £° = in both (3.5) and (4.1)), the time derivative ^2- is also 
solution of (4T) with the same Lyapunov exponent as b°, a = A (= Aq in this case). In this limit we would naturally 



define \E , 2r('?") as ^j-. This property is lost in the more general situation of a non vanishing noise, because £° is not 

time-invariant. What remains true however is the fact that b° and still span the set of "coherent" fluctuations 
which do not affect the shape of the pulse but modify its height and position. 

By turning now our attention to the linear dynamics dual to 
us indeed consider the equation of motion 



l|) we shall come close to the equation (p.6[). Let 



d0 

dr 



'£0, 



(4.4) 



G 



where in order to limit the proliferation of symbols, we keep the same notation 8 for the n ew a uxilia ry fi eld, alth ough 



it is only in particular circumstances, to be clarified below, related to the 8 of equations (3.5) and (3J3). From ( 4.1 
we get 



dO 

d^ 



l M8- d h ({h-h) l / 2 e.B^ 



b° 



The dual dynamics enables one to construct a basis of left eigenvectors ^u(t) (with 1 < i < d) satisfying 
as well as the following orthogonality conditions with the members of the first basis : 

k 

at every time. The vectors ^fu(0) are determined by diagonalizing the adjoint Floquet operator 

1 Ut = exp / l Cdr T+i , 
Jo 



(4.5) 

(4.6) 
(4.7) 

(4.8) 



and enforcing the normalization condition (4.6) at time r = (exp is now an anti-chronologically time ordered 
product). We may note at this point that the first left eigenvector is, in the generic case of non-zero noise, 
the only one to display the scaling behaviour anticipated f or 8 ° according to (3.14), since its Lyapunov exponent 
equals —<j\ = —A. We conclude that the auxiliary equation (3.6) in the restricted manifold of self-similar solutions is 
tantamount to the relation : 



0°(r) =/*i*iiW, 



(4.9) 



where the multiplicative constant fix is nothing but the overlap 8° ■ b° (= /ii^ii(r) ■ \&i. r (r) = /Ui), which wa s sh own 
before to be indeed a conserved quantity. This claim is further confirmed by rewriting the original equation (|3.6[) as 



d8 



1 , 



— - -*M 8-{8B *B0)b - (b • b) d.-[(d h B) l B + B (dh f B)}8 
dr 2 

t\Aa fa d*da\u ("u ul o /a d\ *i 



= -*M 8 ~ (8 ■ B *B8)h - (b • b) 8.(d h B) l B8 . 
Putting back superscripts and reintroducing £° by using ( 3.12| ) , we arrive at 

b° 



(b° ■h°) 1 ' 2 d b {8 Q -B?), 



(4.10) 



which shows that 8° obeys the dual dynamics defined by equation ( iA ) 



Having interpreted (3.6) as a condition of self-consistency for the conjugate momentum 8 expressed by ( |4.9| ), we 
could contemplate the following Newton- like procedure for catching numerically sel f-sim ilar instantons. First make 
a guess for £ in the form of a traveling wave (£ m (r + T) = T+i ^'"(t)), integrate (3.6) forward in time in order to 
determine the asymptotic traveling state reached by b in the prescribed configuration of the noise. Then compute 
i&n from the diagonalization of t Ur (or from running ( |4.4| ) backward in time in order to let emerge the eigendirection 
of lowest growth rate), employ (4^6) for producing a new configuration of 8 (and thereby £) and iterate this loop 
many times at a fixed value of the overlap [i\ until convergence is achieved. However two major difficulties call for an 
improvement of the method : they both have to do with the stability of the trajectory upon time reparametrization. 
First we do not know the speed (or the inverse period T -1 ) of the final traveling wave which must carry together b and 

Therefore, when performing the first step of the iterative loop, we must allow continuous time reparametrization of 
our Ansatz for the nois e in order to fine tune the speeds of the two pulses formed by b and £ and let both terms in the 
right hand side of (3.5) be always relevant. It will be explained in the next Subsection how this goal can be achieved in 
practice. The second difficulty is much more serious than the preceding one and in the wa y of getti ng ar ound it resides 
perhaps the most tricky part of this work. The point is that traveling wave solutions to ( |3.5| ) and ( |3.6[ ) may perfectly 
have a non zero pseudo-energy H, while we are looking for the particular ones with 7i = 0. We shall be able to fulfill 
asymptotically the two conditions h.8 = /ii and Ti = 0, if and only if our iterative guess for 6 is constructed within a 
two-dimensional space rather than a unidimcnsional one as in the naive proposal made above. For this purpose, we 
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are going to embed the linearized dynamics (4.1) into a new one which admits the time translation mode as a true 
eigenstate, of the same growth factor A as b°, restoring thereby the symmetries present in the absence of noise. We 
shall do that in the most economical way, from both formal and numerical points of view, by substituting $2r = 
to 1 4 , 2r, i-e., the eigendirection along which the fluctuations of b around b° are the less stable. 

The left eigenvector \&2( is first rescaled as &21 = .t,^^, — (which makes sense as long as *&2i i &2r 7^ 0, a condition 

always found to be satisfied in practice), so that 3>2Z has a unit overlap with 3>2r, while being orthogonal to all other 
right eigenvectors with i 2 (which will be noted 3>i r from now on). One then considers the modified linear 
dynamics : 



dSb 



CSb = CSb + (b° ■ b ) 1 / 2 (3> 2 ; ■ 8b) B ■ 

dr dr 



It is easily checked that <J> 2 r = ^r- obeys (4.11), since equation (3.5) yields upon time derivation : 



^ = + (b° • b°) 1/2 B d -f=l^r- 

dr dr 



(4.11) 



(4.12) 



It is also trivially seen that the other vectors 4>j r = \l>j r for i 7^ 2 keep the same evolution under ( 4.11 ) as under (|4.l|). 
The dual dynamics reads now : 



f £9 = - l ce - (b° • b ) 1 / 2 (e ■ b ^)*2i 



(4.13) 



It leads to a new family of left eigenvectors &u, dual to the direct basis, whose second member $2; has been defined 
above and the others relate to their original counterparts as 



21 ■ 



(4.14) 



Although the $i;'s were introduced as a rather formal trick, it should be emphasized that 3>i; and $2/ have an 
appealing physical meaning. Parameterizing a perturbed trajectory for b as b = e sln K T ) b°(r + St(t)) + Sbi nc , where 
the "incoherent" part of fluctuations <5bj nc is bound to be a linear superposition of the less dangerous modes 3>j r for 
i > 3, one has indeed, to linear order in (5b, 

S\nb = <f> u ■ 5b , (4.15) 
5T = <f> 2l -5b. (4.16) 

These two relations will be useful in the computation of quadratic fluctuations to be presented in Section They 
show that by projecting out the multidimensional fluctuation field 5b onto the two vectors <&u and one has 
access to the most relevant part of it affecting respectively the amplitude and the time delay of the pulse constituting 
the instanton. In terms of <&i/ and $2/, the self-consistency condition ( |4.S| ) for 9° together with the requirement of 
zero pseudo-energy TL take the following form : 



where 



0°(r) =/ii*i/(r)+/i2(r)*2i(r), 



^(r) = le°.B t B0° = ie.e. 



(4.17) 



(4.18) 



Since from (|4.17p, h 2 {t) = *2r-0°, and from the equat ion o f motion (p|, * 2 r ■ 0° = H + h$° ■ $°, the relation (E3 



is just a way of restating H = 0. That (4.9) implies (4.17) results from the general link between ^fu and G>u (see 
(4.14)). The reverse is true only under the supplementary condition of constant 7i or ^(t) — C te + g£ • £ , which is 



guaranteed by (4.18). It is proven by checking that in that case 9 (r), as given by (4.17), obeys, as it should, (4.4) 



df_ 

dr 



= - f C9° 



dfJ.2 
dr 



2/ 



*L9 a - (b° • b°)^ 2 (9 ■ B ^)<S>2i + ^ * 



dr 



dr 



21 



21 ■ 
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The great advantage of (4.17) and (4.18) with respect to ( |4.9| ) is that this couple of equations lends itself to iterative 
procedures leading inexorably to a fixed point of zero pseudo-energy, a task seemingly out of reach before. There is 
some unavoidable arbitrariness in the construction proposed here, concerning in particular the definition of the vector 
3>2i, about which the reader may feel a little uncomfortable. We suspect that these unwanted features do not affect 
the final results since the original equations to be solved as well as the corresponding conserved quantities all have a 
clear mathematical definition where $1/ and $2/ merge together into ^n. 



B. Practical implementing and results 

The action density so (z) could be computed successfully for the three stochastic models defined in Section || using 
the iterative scheme outlined before. The shell lattice was first mapped onto a circle of d sites, with d typically ranging 
between 20 and 30. Finite size effects turn out to be completely negligible at such lengths of the chain, due to the 
strongly localized structure of the instantons. To start the computation, we make a guess for both the unit field 
^( r ) = (b b/v^r) an< ^ ^ e n °i se £( r ) called henceforth C m (r) and £ m (r). They are such that 

C» +1 (r + T in ) = C»(r), ^ +1 {r + T in ) = C(r), (4.19) 

and 

C™ +d (r) = C*(t), C+d(r) = C(r). (4.20) 

Furthermore, the noise is normalized in such a way that the overlap b • 9 = C • B _1 [C]^, takes on a prescribed value 
111 held as a control parameter during all the steps of the computation. A possible and convenient choice would 
be for instance C m (r) = C°(r), where C°(r) is the deterministic solution of scaling exponent zq, and £ m (r) = 
Hi(b° • b°) 1//2 _B[C°]$y (r), where *?;(t) is the left eigenvector dual to b (r), i.e., in the absence of noise. Fig. |l| 
shows how both vectors 4>y an( l l 00 ^ Hloe at a given time. Their shapes will in fact little evolve as we let the 
scaling exponent z depart from zq. 

In order to allow time reparametrization of the trajectory, we first get an estimate of the instantaneous position of 
the pulse along the shell axis in our trial configuration by computing the following quantity 

n™{T) = Y j n[d]{C™(T)) 2 . (4.21) 

n=0 

The notation n[d] recalls that, due to cyclic boundary conditions, the shell index n is now only defined modulo d and 
that in practice a continuous determination of this integer should be adopted close to the center of the pulse which 



contributes mostly to the right hand side of (4.21). One has by construction n m (r + T m ) = n m (r) + 1 (mod d). 
Having recorded n m (r) and £ m (r) during a whole period T m , we integrate forward in time the nonlinear evolution 
equation for C deduced from ( p. 11 ) by projecting out the longitudinal part of its right hand side : 

^ = {N[C](t) + B[C](r)e n (r')} ± , (4.22) 

Note that the subindex T disappeared from the nonlinear kernel Nr because the Ito-drift term being linear in T does 
not matter in the computation of the action to leading order (we shall see in the next Section how to handle it to next 



to leading order). The most salient feature of (4.22) is that the noise configuration is evaluated in relation not to the 
time r but rather to the actual instantaneous position of the pulse. This means that the time t'(t) is automatically 
delayed or advanced with respect to r, according to the recipe 



- in ( T ')=n{T) or = (n 4n ) _1 [n(T)]. (4.23) 



n 



After integrating (4.22) long enough, a new traveling wave state C ou *(r), £ (r) = £ m { T '( r )) wm usually emerge, of 
period T out possibly different from T m and averaged growth factor 



1 



T+T" 



A out = I (N[C out ] + B[C out ]£ out ) .C out dr. (4.24) 



The vectors *&u(t), ^ , 2z( t ) are then identified as the two eigenvectors of t Uj">nt (defined in (4.8)) of smallest (real 
negative) Lyapunov coefficient and the corresponding «&k(t), &2i(t) constructed as linear combinations of them 
obeying for all times the following relations 
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iih ouf 

* li -b°"* = * 2J .^_ = l 

dr 

dh out 

® 2l .h out = <t> u .—— = 

dr 

Finally, the trial noise configuration is renewed as £ m (r) = (h out ■ b "*) 1 / 2 t B[C out ]9 m with 

r"(r) = /xiShM + M2(r)* 2i (r), (4.25) 
where ^(t) is determined upon imposing the condition of zero pseudo-energy on the trial solution (6 m (r), b out (r)) 

m . N[C out ] + i(b ont .b OM *)6> m • B[C out ] f B[C out ]0 <n = (4.26) 

After setting b m (r) = b ou *(r), we are ready to repeat the operations described above as many times as needed 
until a fixed point of the transformation (such that b out (T) = b m (r) and £ m (r) = £°™ (t)) is reached, which solves 
the problem. The good stability properties of the algorithm, as well as its iterative character, authorize a rather 
unsophisticated handling of issues raised by the time discretization. As required by the Ito- convention, the equation 



of motion (4.22) was integrated using a first-order Euler scheme with a time step At = about 350 times smaller 



than the period. The time r 1 was approximated as the multiple of At making the relation (4.23) best satisfied. 
Similarly no higher order interpolation scheme was devised for estimating with accuracy the output period T out 



it was again simply approximated as the multiple of At making periodicity conditions (4.19) best satisfied for the 
output. However the time step At was changed at each iteration of the loop so to maintain the time resolution N 
constant. The efficiency of the method was greatly improved, when seeking solutions of exponents z far from Zq, by 
increasing z gradually (through an increase of the control parameter /ii) using as first guess solutions of lower but 
close scaling exponent z' . In this way convergence toward satisfactory self-similar solutions (of exponent z varying by 
less than 10 -5 under iteration and pseudo-energy TC < 10 -5 ) was attained in no more than 20 iterations. 

We turn now to the presentation of our results. Fig. |^ shows the action density Sq{z) for the three models (i), (ii) 
and (iii) defined in Section [n| Values of T were adjusted in order to provide the same curvature of sq(z)/T at the 
bottom of the curves, reached evidently at z — zq. We see that the variations of the zeroth-order action get sharper 
on the z > zq side (the only one displayed in Fig. ||), as one goes from model (i) to model (iii). Figures f| || and ||, 
referring respectively to models (i), (ii) and (iii), show the normalized coherent field C and the random force B£ at 
increasing values of z (0.75, 0.85 and 0.95). In all cases the random force is found to be negative at the leading edge 
of the pulse, in agreement with the physical picture advocated in Q : growth can be enhanced only by frustrating 
the energy transfer processes. For model (iii), the coherent field itself gets negative at the forefront : noise in that 
case just helps to prepare the system in an initial condition consisting of a pulse and a negative well in front of it, 
which then collide. An interesting upshot of our computations is that models like (ii) or (iii) involving only a local 
coupling between b and £ escape the disaster met in the framework of model (i), namely a cross-over toward an 
asymptotic linear growth of sq(z) with z, already perceptible in Fig. ^|. Such a behaviour forbids the existence of 
velocity moments at arbitrary orders and is, thus, clearly undesirable. It turns out that the whole shape of so(z) for 
model (i) can be pretty well understood from an adiabatic approximation, which is carried out in the Appendix [A|, 
where solutions of arbitrary scaling exponents are constructed using adequate time reparametrizations and dilations 
of the deterministic solution of scaling exponent zq. The validity of this approximation for the model (i) is somehow 
obvious from Fig. ||, where it can be checked that instantons keep indeed almost the same shape, even for quite 
sizable variations of z. Its failure in models (ii) and (iii) is probably due to too strong deformations of the solutions 
as z increases, again suggested by Figs. ^ and |[ The full non-linear treatment of the problem proposed in this paper 
was however necessary to reach this quite fortunate conclusion. 



V. THE EFFECT OF QUADRATIC FLUCTUATIONS 
A. Formal considerations 

In order to compute the first order (in T) correction to the action per unit cascade step (si(z) in the expression 



(2.8)) of the density of probability P n (z), we have to expand the MSR action up to quadratic order in fluctuations <5b 



around the extremal trajectory b°(T) of scaling exponent z and then sum over them in a way which will be explained 
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below. Since typical fluctuating paths are not differentiable but rather behave as Wiener paths with derivatives of 
the order i, we shall stick to time-discretized expressions in all the following manipulations of the path integral. For 
the sake of clarity, the superscript referring to the extremal trajectory in the previous Section will be taken away, 
whereas b^, 6i will be short-hand notations for b°(T,; = i At), 6°(Ti = i At), where St is the (small) time step used 
in the discretization. 



We start from (3.1) and the representation (3.3) of the effective action S[b, p], expand it to quadratic order in both 
fluctuations Sp and Sb and then integrate out the fluctuations of the auxiliary field. To begin, the time ry during 
which we let the system evolve will be equal to the time nT needed by the ideal instanton to perform n steps along the 
shell axis. The ideal initial and final configurations b in — b and = b^v describe then a pulse centered successively 
around the shells of index and n. To quadratic order in deviations from the extremal trajectory, the probability of 
joining the perturbed endpoints bj„ = bo + Sbg and bf = bjv + Sbjy in the time r/ take the following expression : 



P(b + <5b , 0; b N + Sb N ,T f ) = e~ nso ^/ r J VSb exp (-£S[5b] - S 2 S[Sb}) , 
where the measure of integration is defined as 



VSb 



( 1 



\2hTAt) (b • b ) d / 2 | detBo | JJ 
the linear variation SS reduces to the boundary term 



dSbi 



(bi • b, ; ) d / 2 I dctB 4 



(5.1) 



(5.2) 



ss[Sb] = -{e N ^ ■ Sb N - e ■ Sb } . 



and the quadratic one reads 



S 2 S[Sb] 



At 
2T 



N-l 

E 

i=0 



Sb, 



At 



AiSbi) 



Sbi.ViSbi 



The drift- and potential- terms in eq.(5.4) are found to be given by the following relations 

ti 



A,Sb, = M l Sb l + {Sb, ■ dbj ((bi • bi)Bi %0i 



and 



(5.3) 



(5.4) 



(5.5) 



Sbi.ViSbi = -6i.(Sbi ■ db^MiSbi - -0i.(8bi ■ 9b,) 2 ((b, • b z )B, %6 t ) , 
where M is, as in previous Sections, the Jacobian matrix of the nonlinear kernel Tg^rrj- 



(5.6) 



Our ta sk i s to perform explicitly the integration over fluctuations £bi, • ■ ■ , Sb^-\ at intermediate steps in the path 
integral (5.1). First it will be convenient to get rid of the anisotropy of the "mass" tensor acting in the kinetic term 
of S 2 S |21|] by switching to normalized fluctuating fields defined as 



Sb, = (b.-b,) 1 / 2 ^ 
In this way the measure in the integral transforms into 

<yy. i=jv-i 



in — J —r-z — 2?h with 

(bo-bo) d / 2 [detBol 



Vh 



( 1 



V27rrAT 



(5.7) 



(5.8) 



In performing this change of variables in S 2 S, we must pay attention to the fact that Sbi+i — Sbi, as well as hf+i — h^, 
are potentially of order At 1 / 2 . One finds that, up to 0(At 3 / 2 ) corrections (negligible in the continuum limit), 5 2 S 
becomes : 



S 2 S[h] 



At 
2T 



AT-l r 

E 

i=0 



jl/2 ^hj+i 



An 



b 1+Vt b 1±1 +b 1 ^ 



(5.9) 
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where 



Qi = 



bi ■ bi 



1/2 



H+l ■ L»t+1 

b, • b; 



1/2 



Bj 1 i?i+i 



(5.10) 



-A^ — B- j\.<i Bi , 



(5.11) 



and 



At 



V; = {bi-b^BiViBi, 



(b i+1 • bi+i) 1 / 2 ^! + (bi ■ b0 1/2 sj " Ubi+i.bi+i^Bi+x - (bi ■ b, ; ) 1/2 B 



A few simplifications are now in order. On one hand we expand Qi as Qi — 1 + SQi + 0(At 2 ), where 



Pf+i • bj+i 
b, ■ bi 



1/2 t 



S i+i '-^i 1 + B i 1 Bi+i 



(5.12) 



(5.13) 



(5.14) 



is of order At. Tracking O(At) terms in S 2 S, we may replace Qi by 1 everywhere in (5.9), except for the kinetic 
term -p (h i+ i-hi).Q^{h i+1 -hi) ^ w j 1 j c j 1 accor di n g to standard computation rules JTs| in path integrals can be reduced to 

W < ' h ' + AT h ' ' > — ^ TrJQj. On the other hand, exact ways of tracing out quadratic forms like 5 2 S are available only once 
they are written in terms of the mid-point field h ' + ^ +h ' rather than (which amounts going back at this stage to 



the Stratonovich prescription). Again to order At, the substitution of h ' + ^ +h ' to h^ in (5.9) is harmless except in the 
cross-product term — jr(h i+ i — hj) ■ A'^hi which gets into — jr(h i+ i — h,-) ■ A[ — " H — ^~ Tr.4-. 

Putting things together, we arrive at the following expression for the transition probability in the neighborhood of 
the instanton : 



P(6b m -> Sbf,T f ) =| det^(r/) | e -^ r e" 1 ^ e^«-^-0o-<5b,„ } ^ ^ 



where we defined 



N-l 



J ^ T f) = \ E Tr ^ + ^ E Tr ^ - dAr / ' 

i=0 

and the reduced path integral Z[hi n — » hy, t/] as 

Z[h in -► h/,r/] = 



AT-l 



i=0 



(5.15) 



(5.16) 



hjv=hf 



X>h exp-(5 2 5[h] 



h =h iB 



with 



AT-l 



i=0 



At 



hj+i + bi hj + i + hj 



and 



(5.17) 



(5.18) 



(5.19) 



One easily shows that in the At — > limit 

^ Tr<5Q 4 -> Tr 



In 



(bjy ■ bjv) 1 ^^ 



A* 



(b -b )l/2 Bo 



dAr f , 
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since the periodicity of the instanton implies that det-Bjv = det-E?o, while 



ArJ> A - [ {^.[{b ■ b)B *B6] ^^}dr. 



However the ensuing expression for Zi(t/) is not yet complete to order O(T ). This is because we discarded the 
O(r) Ito drift-term in our computation of extremal trajectories. The small deviations induced by this term may be 
neglected in Xi(r/) and — lnZ[hi„ — > hj,rj] which are already first order corrections in a T-expansion but they must 

be taken care of in the zeroth order term. Rather than n ° _f , it should read Sr[b r ], where the index T denotes 



a quantity or field evaluated in the presence of the Ito-term specified in (2.1C). To first order in T, we can take 
advantage of the extremum property of bp and write down 

SrK] - S[b°] = SrK] - S r [b } + S r [b } - S[b°] 

~ 9 fh°l 9fh°l - 1 r 0° N ( b °]- N r[b°] f]T 
~S r {b }-S[b ]_- ^ 9 . (b0 ^ bQ)1/2 dr, 



where the last identity just comes from the expression ( p.4[ ) of the action. Rearranging things under the assumption 
(satisfied by each of the particular models that we considered) that the entry Bjk of the matrix B involves homogeneous 
monomials of degree I built up only from components C m with m ^ j (one has I — for model (i) and / = 2 for 
models (ii) and (iii) defined in Section |l|), one finds that Ii(t/) in ( |5.15 ) should finally be understood as : 



MT f )= { ^-^Ar f + ^LJl r^.BWdr. (5.20) 



n 



This preliminary work being done, the discussion will now concentrate on the reduced path integral Z[hi n — > hy, Tf]. 
We restrict first our attention to the case of fixed endpoints h in = hf =0. The instantons found in the previous 
Sections are physically satisfactory only if the quadratic functional 5 2 S'[h] is positive for all {h^} such that ho = hjy = 
(we shall see later on that this is not a sufficient condition in the present problem) . According to standard results of 
functional analysis |l9(| , positiveness of i5 2 5[h] is tantamount to the absence of points conjugate to the origin during 
the whole time interval [0, r/]. Recall that the definition of conjugate points goes as follows : d being the dimension 

of the space (here the number of shells) we construct d initial conditions (h Q Q ' ) , ) such that 

4?=0> h w = At<5 q/3 + 0(At 2 ), 



where j3 is a shell index running as a between 1 and d and let them evolve under the Euler-Lagrange equations derived 
from <5 2 5[h]. The time Ti — iAr is said to be conjugate to the origin if the system formed by the d vectors gets 
degenerate there. One can build a matrix Ui such that — h?, in terms of which the initial conditions read 

U = 0, Ui = At + 0(Ar 2 ) , (5.21) 

while Ui+i (for 1 < i < N — 1) is obtained from Ui-\ and Ui through a matrix Euler-Lagrange equation, derived from 



(5.18) and conveniently cast into the following form 



where 



Pi = - Bi U -±±^± , (5.23) 

can be seen as a matrix momentum. In the absence of conjugate points, det Ui never vanishes except at the origin 
and one gets the following simple expression for the reduced path integral, provided the Ar — » limit is ultimately 
taken : 



d/2 

\2ttt) y/detU(T f ) 



Z[h m =0^h / = 0,r / ]= f-^-V 1 (.-,.24) 
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Details on this result, which may be found in many textbooks on path integrals pO|pl|] , are provided in the Appendix 
[b|. A very nice feature of the proof of the connection of the positiveness of 5^S[h] with the absence of conjugate 
points is that it provides also an efficient way for computing the reduced path integral with h/ arbitrary (but still 
small naturally). Indeed the main idea consists in adding to J 2 S'[h] a boundary term of the form 



1 



N-l 



where W is a symmetric matrix at all times and then selecting the right W such that 5 2 S + S 2 S' becomes a perfect 
square. In order to achieve this task, Wj must be a solution of a matrix Ricatti equation (see again the Appendix 
B for details) which, upon the substitution of a new unknown matrix f/j defined implicitely by the relation (for 
< i < N ~ 1) : 



1 



{WiUi+Wi+xUi+x) 



U i+1 - Ui 



U, + u. 



(5.25) 



is found to be nothing but the matrix Euler equation ( 5.22 ). It follows that the result ( 5.24 ) may be extended to the 
case of an arbitrary final configuration hf as 



Z[K 



d/2 



(5.26) 



where Uh-f) and W(tj) are to be computed by letting both matrices evolve from r = to r = Tf according to ( 5.22 ) 
and (5.25). There are some subtleties about the choice of initial conditions and proper counting of the number of 
unknowns whose discussion we prefer to relegate in the Appendix |b| 



B. A physical definition of si(z) 

We are now in a good position to compute the next to leading order term Si(z) in the expansion of 
— lim„^ +00 ilogP„(z) in powers of T. We shall obtain an estimate for P n (z) by summing over all the trajecto- 
ries which lead to the same growth of the pulse as the ideal instanton b°(r) after n cascade steps. In the small T 
limit, all statistically relevant trajectories remain close to b (r) and we may define unambiguously their "arrival" 
time at the shell of index n as the first time r„ such that 

b(r n ) = b°(T° = nT)+Sb', (5.27) 

where <5b' reduces to a linear combination of stable "irrelevant" modes 4v(r„) for i > 3. Up to multiplicative factors 
growing at most algebraically with n, we can then write P n (z) as the following integral over the arrival time and the 
position of the endpoint 

P n (z) « J P(b°(0) b°(T°) + «5b',r n )5(* lz (r°). < 5b , )<5(* 2Z (T°).(Jb')dr n d5b'. (5.28) 



The density of probability in the right hand side of this expression is known from (5.15) and ( [5.26 ), where Tf should 



be taken equal to r n and dhf equal to £b = b(r n ) — b°(r„). Calling 5t — r n — r° the time delay, we get the following 
relation between £b and ^b', valid up to 0(5t 3 ) corrections : 



^b = SW + b°(r n °) - b°(r„) = 5b' - Sr^(r n ) + ^r 2 ^(r„ ) . (5.29) 



Note that St scales typically like VT 7 in the semi-classical limit. We may thus, to leading order, replace T n by 
r° and identify Sh with Sh' — Jr^ r (r°) everywhere in the integrand of the right hand side of ( 5.2§| ), except in 



y(So(tu) + 9 (r n ).5b) (the combination appearing in the exponential prefactor of (5.15)) which deserves more care. 
It follows from the definition of the action that (again up to 0(St 3 ) corrections) 

S (T n ) = S (t°) + St - \ST 2 ^ (r„°).^(r°) , (5.30) 
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and from fl5.29[ ) (together with the condition u (r°) ■ 5b' = 0) that : 

0°(r») ■ 5b = 5r ^(r°) • 5b' - 5r 0°(r„) • ^-(r n ) + l S r 2 e (r n ).^(r° n ) , (5.31) 

When summing these two equations, linear terms in St disappear as expected and we are left after some rearrangements 
with the remainder, of quadratic order in fluctuations : 

^{S (r n )+e Q (r n ).Sb} = + i ( 5r 2 ^-(r n °) • ^(r") + St ^(r°) • 5b, (5.32) 
db° 

where we set 5b = 5b' — St (t2), so that the time delay St is simply expressed in terms of 5b as St = — 3>2;(T?)-5b. 

The end of our theoretical considerations is reached with the following expression for P n {z) : 

P n {z)^e- nso ^l v , = / — f — — 5(*i,(r°)-h)dh, (5.33) 

V ; y/dctU(T°) J (2^r) rf / 2 V m n > ' V ' 

where #i,(r) = (bO-b ) 1 / 2 *B*i/(t) (so that *' ir h = #i;-5b), W> = T^ + AW, with W defined in the last subsection 
and 

h.AW h = -2(*^ • h) (b°.b°) 1/2 ^r • £h + • h) 2 ^- ■ ^- . (5.34) 

We see that the condition of positiveness of the matrix U must be supplemented by the condition of positiveness of 
the restriction W\ of W to the (d — l)-dimensional space orthogonal to & u , in order to make the instantons found in 
the preceding Section physically meaningful. Provided these two requirements are met, si(z) is obtained as 

si(z)= lim -{X 1 (nT) + ilndetC/(nT) + ilndetW / i(nT)}, (5.35) 

n^+oo n 2 2 

which is the main result of this Section. 

Tracing back all the steps leading to ( 5.35| ), one could object to our starting point ( 5.28 ) the fact that fluctuations 
of the initial endpoint are not taken into account. This could be done at the expense of rather more cumbersome 
formulae for the reduced path integral Z[h(0) — * h(r)] when both h(0) and h(r) do not vanish. However, we believe 
on physical grounds that the blowing-up associated with the instanton washes out any influence of the fluctu ation s 
at large scales on the part of the action scaling linearly with the number of steps n. Therefore the expression ( 5.35| ) 
should be exact. 



C. Practical implementing and results 

The most difficult part of the computation of si(z) lies in the evaluation of detJ7(r) and detWi(r) (as defined 
in the preceding subsection), which requires a good control of all the eigenvalues of these two matrices. Numerical 
instabilities could be avoided for a time long enough to get a p recis e estimate of Tr-ln det UW\ by using the exact 



expression of the matrix Euler equation (5.22) and the relation (5.25) between W and U. In this problem, there are 
two Goldstone modes associated with uniform rescaling and time translation of the instanton : as a consequence, the 
matrix (b° • b°) 1 / 2 i3?7(r) (resp. (b° • b ) -1 t B^ 1 WB~ 1 ) is expected to have two eigenvalues scaling like r (resp. like 
1/t) (in order to obtain the simplest transcription of these symmetries, one has to go back to the original fluctuation 
field 5b = (b° • b°) 1 / 2 _Bh and the change of variable influences eigenvalues of U and W). When restricting W to 
the (d — l)-dimensional space orthogonal to 3>' i; (= (b° • b°) 1 ^ 2 B^u), one looses one of the eigenvalues scaling like 



1/t, so that there remains an algebraic factor y/r in the product \J det U[r) y W% (r) , which we had to take away 
by hand in order to make more conspicuous the leading exponential growth of this quantity. To give an idea of the 



accuracy of our procedure and confirm the soundness of the intricate formula ( 5.35| ) that was proposed for si(z), 
we show in Fig. o the behaviour of various relevant quantities in the case of model (ii) , and for a moderate scaling 
exponent z = 0.8. It is observed in the picture on the top that the logarithmic derivatives of det U and det W exhibit a 
linear behaviour with almost opposite slopes. It can be shown by considering simpler and exactly solvable models for 
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quadratic fluctuations without inter-shell couplings, that this strange effect mostly reflects the stiff variations suffered 
by the noise variance on any shell, as the latter moves back from the leading edge of the instanton to its rear end. The 
argument is presented in the Appendix ^|, since it may help the reader to get a feeling for the order of magnitudes at 
play in both U and W matrices. 

It could be tempting at this level to approximate si(z) as 

lira -^-{Ti + \ detlnW}, (5.36) 
t-^oo dr 2 

This expression would come out if, without great physical justification, it were decided to sum over all positions of 
the final endpoint in the path integral formulation in order to estimate the volume of the basin of attraction of the 
instanton. We found however that the matrix momentum WU develops invariably a negative eigenvalue after some 
time, so that the projection onto the restricted phase space introduced in the previous subsection is a necessary step 
for restoring the statistical stability of the instanton. Further, even before this instability occurs, there was found 



to be a residual t 2 term in IndetWU which forbids any reliable estimate for s\(z) to be deduced from (5.36). The 
picture at the bottom of Fig. |6| shows by contrast that the more precise quantity det U det W\ quickly settles to a 
perfect exponential growth, once algebraic transients have been factored out. Note that the positiveness of both U 
and W\ could be checked in any instance. 

In all the models that we investigated, we found that the first order correction Si(z) to the action takes an approximate 
parabolic shape of positive or negative concavity, centered around a value of z different from zq. In the case of the 
model (ii), the concavity of si(z) is opposite to the one of sq(z) and the maximum of si(z) is reached at a scaling 
exponent z\ w 0.6, significantly lower than zq = 0.72. This means that the minimum of the total action density 
s(z) — sa(z)/T + s\(z) (for values of T such that s(z) remains concave as it should) is displaced toward the side 
of larger exponents, just as a result of fluctuations. The trend is just the opposite for the model (iii), where s\(z) 
presents this time the same concavity as so(z) and a minimum on the left side of zq. Fig. shows the graph of 
sq(z)/T + si(z) that is obtained for the model (ii) and for the value T — 0.58, which we believe to be of some relevance 



for the GOY model (see Section VI) 



VI. COMPARISON WITH NUMERICAL DATA ON THE STATISTICS OF COHERENT EVENTS IN 

THE GOY MODEL 

The systems we have analyzed here were introduced to describe inertial singular structures of the GOY model. To 
check their physical relevance, we first have to define a class of events observed in full simulations of the GOY model 
which are likely to be the best candidates for a description in terms of instantons. It is clear that relative maxima of 
the instantaneous energy flux e„(i) (with e n = £;~ 2 Re{(l — e)Q 2 b n ^2b n -ib n + b n -\b n+ ib n }) are useful observables for 
tracking the passage of coherent structures across the whole inertial range. But their total number is found to grow 

2/3 

with n as k n , due to the acceleration of time scales typical of the Kolmogorov energy cascade. One may consider that 
they develop on tree- like patterns in the (n,t)-plane, which are renewed at each turn-over of the large scales (totls). 
We say that such trees provide a realization of the propagation of a coherent event from shell no to shell n > uq, 
whenever the n — no nodes of the tree closest in time to their supposedly common ancestor on the shell of lower index 
no appear in the order of increasing shell index. In order to discard too weak, and therefore irrelevant, events, we 
imposed that e„ be greater than half the mean energy flux. Figure || shows the logarithm of the histogram of the 
logarithmic amplitudes A n — In le^l 1 / 3 for all relative maxima on one hand and for the restricted class of coherent 
events defined above on the other hand, with hq — 5 (far enough from the forcing range) and n = 11 (well in the 
inertial range). The Reynolds number of the simulation is Re = 10 8 . Statistics have been run over 6 x 10 4 totls, and 
in average there are three "coherent lines" for two totls. We note that the statistics of coherent events is very close 
to log-normal for |e„| > 0(1). 

The effective exponent z of a coherent burst after n — no cascade steps is obtained via the relation A n = A no + (n — 
no)(z — 2/3) In 2. If anomalous scaling is preserved in the Re — > oo limit, the pdf of scaling exponents z should behave 
at large cascade lengths as P n {z) ~ e ~ ns ( z ) j where generically s(z) will present a quadratic minimum at some z*, with 
an expansion around z* that we write as s(z) = a(z — z-A 2 + • • •. The histogram of the variable A n — A no should 
consequently evolve, as n — uq increases, towards a Gaussian shape, whose center D n and variance E 2 grow linearly 
with n and relate to z* and a as 

D„~nln2(2*-2/3), £ 2 ~^(ln2) 2 . (6.1) 

la 

The actual behaviours of D n and — 2£ 2 , obtained from a very high Reynolds number simulation (Re — 10 9 , which 
sets the dissipative scale at the shell index = | ln 2 Re = 23) are plotted in Fig. ||. Error bars were estimated by 
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varying the range of the quadratic fit to the logarithm of the histogram of A n — A no , as well as the domain of initial 
amplitudes A ng used to construct this histogram. It appears that the two graphs are rather far from simple straight 
lines : this is especially true for the variance, whose graph from concave gets convex beyond the shell index 15. From 
the investigation of lower Reynolds numbers Re — 10 8 and Re = 10 7 we could deduce that this transition occurs 
at a shell index n c always of the order of rid — 8 and defines a clear-cut boundary between the inertial range and a 
surprisingly wide pre-viscous range. We believe that the direct action of viscosity on intense bursts starts to show 
up only at the shell index rid — 3, beyond which the local slope of the -D„-graph ceases to vary and points to a value 
of the average growth exponent of coherent structures precisely equal to zq. The fact that n c lies rather far from rid 
means that the cut-off imposed by viscosity exerts a long-range influence on the statistics of the random environment 
seen by a coherent structure. A decent linear regime for both the drift and the variance is observed in the range 
15 < n < 21, from which we get the two estimates a — 29 ±4 and — 0.74 ± 3. 10~ 3 . Assuming that the fitting range 
10 < n < 16 provides the best clue to the asymptotic scaling of the inertial range, one gets the second set of values 
a = 45 ± 6 and z+ = 0.75 ± 3.10 -3 . It appears that the physics of the previscous range is quite well reproduced within 
our modelling (ii) of the incoherent background. By choosing V = 0.58 (a value hopefully small enough to fall within 
the range of validity of semiclassical approximations), one obtains, as Fig. |?] shows, an almost perfect parabolic shape 
of s(z) = —(Sq(z)/T + S\(z)), with a maximum reached at z* = 0.74 and a curvature a = — 2s"(z±) = 29. However 
the parameter T cannot be adjusted so as to account for the higher values of a and z* characterizing the inertial 
range. It would seem that in this range of scales it gets necessary to assume some bias in the incoherent fluctuations 
boosting the increase of the renormalized value of z*, while keeping the noise width small. 



VII. CONCLUSION 



We have developed a general scheme for computing numerically self-similar instantons in scale invariant stochastic 
dynamical systems. As concerns the physics of the GOY model, we believe that the bunch of results presented here 
give a strong support to the relevance of an approach focusing from the outset on structures in order to understand 
intermittency and treating the rest of the flow as a noise of weak amplitude. In particular the trend toward log-normal 
statistics of coherent structures is nicely recovered. The detailed study of various stochastic extensions of the GOY 
model shows that the resulting pdf of scaling exponents of singular strucures is very sensitive to the hypothesis made 
on the coupling of noise to the velocity gradient fiels. 

We hope that our approach will be useful for attacking the 3D-Navier-Stokes dynamics along similar lines, once an 
adequate decomposition of the flow into coherent and incoherent parts will have been introduced. An application of 
the method to the Kraignan's model of passive scalar advection formulated on a lattice of shells has already been 
attempted p2| . It has given encouraging results with regard to the validity of a semi-classical analysis, even in 
situations where a small parameter (like Y in the present problem) is missing. 
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APPENDIX A: 



We carry out in this Appendix the adiabatic approximation alluded to in Section IV B. We look for self-similar 
instantons within the restricted manifold of configurations of the type 

b(r) = e*^b°(f ), (Al) 

where b (t) is the deterministic solution of scaling exponent zq, f can be thought of as a "proper" time referring 
to the actual position of the pulse. The two variables r(f) and x(f) parameterize then local cha nges of speed and 
amplitudes of the pulse, which keeps the same shape as in the absence of noise. Note that if (Al) is to represent a 
self-similar instanton of scaling exponent z ^ zq, x(t) must obey the constraint 

x(f + T )-x(f) = (z-z )logQ. (A2) 

We plug now the Ansatz ( |Al[ ) in the equation of motion (2.11) and project it onto the two directions $i r = b° and 

*2r = ^r- We get by doing so 



dfdx 

CLT OT 

^-1 = * 2I ..B£. 
dr 



(A3) 
(A4) 



If the other dimensions of the configuration space are neglected, (A3) and (A4) form a closed two-dimensional system, 
which may be rewritten as 



dx 

Tr =Cl ' 
dr _ 



where correlation functions of £1 and ( 2 read 



dr 



(C<(r)G(f )> = Kr ~r) 1 < < 2, 



with 



3 1 



The Gaussian action density associated to one cascade step within this restricted stochastic system is given by 

r T 



1 C dr 

2 1 c<y-M> 



(A5) 
(A6) 

(A7) 
(A8) 

(A9) 



Once expressed in terms of the diffusing variables x(f) and r(f),it becomes 

s = \ £ dfii^r 1 [i 2 ^ 1 )!! + 2x(V-') 12 + [V-% 2 ] + ^(V- 1 ) 

-[xiV-^iz + iV- 1 ^]}. 
The extremization of s with respect to 4? leads to 
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J = [(i 2 ^ 1 )!! + 2x + (V- 1 : 



22 ) V22 



I 1/2 



(A10) 
(AH) 



22 



1 1/2 



and to an effective action for the remaining variable x 

Seff[x(r)} = [ T ° df{[(x 2 (V- 1 ) 11 +2x(V- 1 ) 12 + (V- 1 ) 22 )(V- 1 ) 
Jo 

-[xiV-^ + iV- 1 )^]}. (A12) 

Assuming the coefficients Vu, Vi 2 and V 22 to be almost constant inside the time interval To, one deduces an analytic 
expression for sq(z) from s e // by just replacing in the integral x by (z — zq) logQ/To (which follows from QA2|)). One 
gets in particular for large enough z — Zq, 

s (z) ~ (z - z )\ogQ {((y- 1 ) 11 (V- 1 ) 22 ) 1 / 2 - (V- 1 )^} , (A13) 

i.e, a linear behaviour as observed for the true solution of model (i). 
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APPENDIX B: 



We derive the formal expression of the reduced path integral Z[0 — ► h/,r/] given in eq. (5.26) of the Section [y] of 
the text. The proof is presented in many textbooks on path integrals but, as far as we know, always using a continuous 
definition of time. This leaves some ambiguity in the right equations that matrices U and W should obey, once time 
is discretized for computing purposes. We found that this issue is crucial mostly for evaluating W and preserving its 
symmetry properties. This is why we feel it useful to show how every step of the proof given in the continuum limit 
receives an exact transcription in the discrete time case. 
We start from the quadratic functional : 



5 2 S[h] 



At 
2T 



i=N-l 

E 

i=0 



H+l 



h J+ i + hi 



H+l 



h, . ., hj + i 



V 

■ " 1 



(Bl) 



Upon the addition of the boundary term — ^ J2i=o {hj+i-Wi+ihi+i — hj.Wihi}, it becomes without any approxima- 
tion : 



i=N-l r 

i=0 
h i+l 



Q 



1/2 



An 



H+l 



(B2) 



where 



and 



V, - V' + t B i B i 



»+i 



Wi 



Wi + W, 



Wi + w l+1 



(B3) 



(B4) 



We see that <$ 2 <S[h] reduces to the time integral of a single (positive) square, if and only if Wi is such that for all 



< i < N- 1 



Vl+'BiBi 



W, 



i+l 



w, + w t 



i+l 



)Q7\Bi 



Wi + W i+1 , 



0. 



(B5) 



We note that (B5) forces Wi to remain symmetric for all time, provided that Wq (arbitrary at this stage) is chosen 
to be such. In order to solve the above Ricatti matrix equation, one makes the following change of matrix variable 



(Bi 



W t 



c/m } = Qi u i+1 



Ui 



At 



(B6) 



^From the expression (B3) of Qi, it is easily shown that (B6) is equivalent to the equation ( 5.2 5| ) quoted in the text. 



Furthermore, by multiplying both sides of (B5) by 



Ui + Ui+i 



on the right, one gets for < i < N — 1 



Wi+jU^-WjUj 
At 



= — B, 



U, 



i+l 



U, 



At 



(%Bi + Vl) 



Ui + Ui+x 



(B7) 



By half-summing the two relations yielded by (B7) at subsequent values i — 1 and i of the temporal index (with then 
1 < i < N - 1), and using Q after noticing that w^u^-w^u^ = (W t+1 u, +1 +w t u^-(mu,+w^u^) 
eliminate W and check that U obeys the matrix Euler equation Q5.22 ), as promised in the text. 



one can 



So far, we have proven that, as long as the matrix U may be inverted, the positiveness of 5 2 S is guaranteed, since in 
that case the matrix W exists at all times (from (B6)) and allows one to transform the initial quadratic form into the 
time integral of a si ngle squar e. W e show now how these matrices lead to a compact expression of Z[h.Q = — > h^, r/]. 
We first note that (jB6j) and (|B7|) provide 2N relations for 2(N + 1) unknowns {U Q , ■ ■ ■ , U N }, {W 0l ■ • ■ , W N }. This 
gives much freedom in the choice of Wq and Uo. In the particular case of ho = 0, it is convenient to set Uo = 
and WqUq = 1 (which should be understood as the limit as e — > + of Uq = e and Wq — e _1 , so that Wq is indeed 
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symmetric). A quick inspection of flB(j ) and (B7) reveals that U and Wi behave then respectively as iAr and (iAr) 1 
to leading order in At for 1 < i < N . One has for instance the exact result W\ = (At) -1 - Bo+ * Bo + ^BqBq + V ). 
Recall that the quantity we wish to estimate reads now : 



Z[0 -> hf,r f ] = e -^ h f WNh f J VhS d (h N - hfje-^^o^'+i-Q'^'+i 
where we defined the new field i/> i+1 (for < i < N — 1) as 



= h i+ i - hi - ArQ i (Bi+ ) 



and the measure of integration T>h as 



N-l 



ph^n 



dh i+1 



l { (2nTAT) d / 2 



i=0 



With the help of (B6), the transformation (B9) may be rewritten as 

, u u fUi+i-UA (Ui + Ui+i y , 

V'i+i = h s :+i - hi - 



hi + hi + i) . 



(B8) 



(B9) 



(BIO) 



(Bll) 



This relation is easily inverted by setting h; = Uid, which gives 

' Ui + U i+X U l+ i -Ui ( u + U i+X \ _ x U i+ i - u 



Ui+i 



U + U+i 



2 

-1 



(Ci+i - Ci) 



(Ci+l Ci) i 



so that we deduce (under the assumption ho = 0), for < i < N — 1, 



hi+i — Ui+i 



E 

3=0 



Since hi+i is linearly related to the ^ - +1 's of index j lower than i, only the diagonal blocks E/j+i 
the Jacobian of the transformation and one has : 

5h l+1 ^ det(Z7i + Ui+i) 



JY 



n 



det 2C/f 



(B12) 



enter 



(B13) 



To enforce the boundary condition h^v = h/ at time Tf in terms of the new variables we introduce the usual 
integral representation of the 5-function : 



5 d (h f -h N ) = J 



da 

(2tt) c 



i _ia.[h / -E, 2 l- I ^ ±i V' i+1 ] 



After performing the Gaussian integration over the i/'j's, one arrives at 

N-l 



z[o^h f ,T f ]= n 



4=0 



det(t7j + U i+ i) 
det 2U 



1 



e -^h f .w N h f x 



det Qi 



da 

(2tt)° 



iQ:h /e --^jEa.GQ: 



(B14) 



(B15) 



where 



G = Un 



N-l 

E 

2=0 



(B16) 
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This awkward non local operator G is greatly simplified when the singular initial condition already mentioned : 
Uq = e, Wo = e -1 with e — > + is adopted. In that case, G is completely dominated by the first term of the series in 
the r.h.s. of (B16) which diverges as e -1 and, to leading order in e, one has 



1 



Gf~_ U N {U^W - U U 



— iw — l*rr— 1\ 'jy^ 



(B17) 



The summation over a can then be done and, since G 1 vanishes in the e — ■» + limit, one gets 



ZfO-h/,7-,] 



/ 1 



V 27rrAr 



det(t/ + t/i) 
det Un 



n 



d et(t/j + ?7 i+ i 
det 2J7i 



det Qi 



(B18) 



Note that all the manipulations presented in this Appendix were devoid of any approximation. It is finally a 
straightforward matter (details will be skipped here), to check that in the At — > limit, the infinite product in front 

v d/2 ! 



of the exponential in (B18) reduces to ( 27r rAr ) Vdct u 
in the text. 



, making thereby (B18) identical to the result (5.26) quoted 



APPENDIX C: 

We consider the following quadratic action : 

= f dr^f- (CI) 

where, in order to mimic the stochastic models studied in this paper, the variance a n (r) evolves on each shell as 

a n (r) = (b -b )(C n _ 2 C° n _ 1 ) 2 , (C2) 

i.e., as the variance of noise, (b° • h°)B nn , for the model (ii). In (|C2|), b°(r) may be thought of as the deterministic 
self-similar solution, and a n (r) can therefore be cast into the form 

a n (r) =e 2A ° T ~a(T~nT a ), (C3) 

where the function a(r) satisfies 

a(r + dT a ) = o(t), (C4) 

because of the periodicity of the shell lattice. Let us assume that the instanton is centered around the shell of index 
n = at time r = 0. At its leading edge (n > 0), C n decreases very abruptly as exp— cr™, with r=(\/5— l)/2 and c 
a constant of order 1. This essential singularity comes from the necessity of balancing with the dominant term 
Q 2 (l — i)C n -iC n -\ of the non-linear kernel of the GOY model in this range of scales. In the trail of the instanton 
(n < 0), one has a much smoother behaviour C n ~ Q nz ° = e nA ° T ° . When the shell lattice is periodized, the leading 
edge and the tail of the instanton have to be glued together and the locus of matching, as well as the residual amplitude 
of C n at that place, will be imposed by the side supporting the slowest variations of C n . We conclude that in a cyclic 
chain containing d shells, most of them reside in the exponential tail of the instanton, so that we may write (again 
under the hypothesis of an instanton initially centered around the origin n = and with a shell index n defined 
between and d — 1) 

a„(0) ~ e 4 (™-<*M° r °, (C5) 

Thus, the range of values spanned by the function a is very large and scales with the total number of shells like e 4dA ° T ° . 
For shells on the exponential ramp, a„(r) first decreases exponentially in time like eT 2A ° T (because C n decreases like 
e ° T in this region) and goes by a sharp maximum of order e 2nA ° T ° at the time t„ ~ raXb when the center of the 
instanton reaches the corresponding shell. Then it starts again to decrease exponentially. 

Having understood these basic dynamical features, we can compute the matrices U and W introduced in subsection 
VA for the quadratic action given by (|Cl|). To make contact with the normalized field h used in the real problem, 



we switch from the variable x n to the variable y n — x n 

I \[&~n- This transforms the original action into 
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d-l 



S2[Vn] = \Y. / dr(y„ + i%„) 2 . (C6) 



2 n=0 ^° 2 °" 

Since there is no inter-shell coupling, the matrices U and W are diagonal in the shell index. The extremization of 
S[y n ] with respect to y n leads to the couple of first-order differential equations 

Pn = Vn + -~ Vn, (C7) 

Pn = \ —Pn- (C8) 
2 a n 



One has simply U nn — y n and W nn U nn = p n . The solution of Eqs. flC7| ) and ( p8| ) under the initial conditions 
y n (0) = and p n (0) = 1 is 



/ \ -y/an(-r) , . 

Pn{T) = = , (C9j 

v a «(°) 

y n (r) = -jl=^l= f a n {r')dr' . (CfO) 



As long as r < t„ = nTo, the instanton has not passed through the shell of index n and the integral in the right hand 
side of (|C10| ) is dominated by the neighborhood of the lower bound r = 0. We deduce that for indices n > t/Tq 

Unn(T) - (Cff) 

W nn (r) =2A a e~ 2AoT . (C12) 

By contrast, when r gets larger than r„ (by some units of time To), the integral in the right hand side of ( |C10|) is 
dominated by the neighborhood of the time r n where a n takes its maximal value. We get 

2m(t) 

a n (0) 

where / is a number of order 1. It follows that for indices n < t/Tq, 

U nn {r) - /e 2(d -" )AoTo e AoT , (CI3) 
W m (r) - r 1 e 2nAoTo e~ 2Ao \ (CI4) 




where we used the estimate (C5) for a„(0). Since W nn U nn ~ e~ A ° T for n > t/T and ~ e 2dA T 0e -A T f or n < T /y 0; we 
conclude that det Wi7 increases exponentially like e dj4or . But this property is not shared by det U or det W c onsi dered 
individually. Indeed, since the number of shells crossed by the instanton increases linearly in time, eqs. (Cll) and 
( |C13T) show that 



I 



r/T 



del r -(^-) e 3AodT e" Ao %, fC'LVi 



while from eqs. QClg ) and ( |C14| ), 



-r/To t2 



This argument captures apparently a good part of the physics of fluctuations around a moving self-similar system, 
though badly treating hybridization effects between neighbouring shells. It also explains how large (resp. small) 
numbers are generated in the spectrum of the matrix U (resp. W) and why in practice one does not have much 
freedom in the choice of the total number of shells d. 
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FIG. 1. On the upper picture (resp. lower) we plot the configurations at a given instant of the right (resp. left) eigenmodes 
in the subspace of maximum Lyapunov exponent Ao around the deterministic solution. 
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FIG. 2. Evolution of the normalized action per unit cascade step so(z)/V , as a function of the effective scaling exponent z, 
for the three models studied in this paper. As a guide for the eyes we show the parabola (solid line) "tangent" to the curves 
at the deterministic minimum zo = 0.72. 
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FIG. 3. Configurations of the normalized coherent field C and the random force B£ for three instantons of exponents equal 
to z = 0.75, 0.85, 0.95, obtained with the model (i) (according to the nomenclature defined in the text). Note that the C field 
is only slightly deformed as z increases. 



2G 




FIG. 4. Same as in Fig.g but for the model (ii). Note the more pronounced deformation of C upon increasing z, with respect 
to the previous case. 
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FIG. 5. Same as in Fig.^ but for the model (iii) and only instantons of exponent z = 0.75, 0.85. The physical field now 
changes its sign at the leading edge of the pulse. 
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FIG. 6. Test for the convergence of several relevant quantities entering the calculation of the function Si(z). We show the 
case of the model (ii) at z — 0.8. 
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FIG. 7. Graph (stars) of so(z)/T + si(z) for T = 0.58 and model (ii). The parabolic fit S qw d = a ( z — z *) 2 yields = 0.74 
and a = 30. 
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FIG. 8. Histograms of the energy flux in log-log plot, involving all relative maxima of e n or only those associated with 
coherent events. The shell index n = 11, the Reynolds number Re = 10 8 and the number of totls is 6 x 10 4 . 
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FIG. 9. The two quantities D n and — 2E^ (encoding the Gaussian central part of the histogram of the growth variable 
A n — A no ) vs n. The Reynolds number Re — 10 9 and the dissipative shell has index rid = 23. The two pieces of straight line 
show the linear fits that were used to extract the values of z* and a in the pre-viscous range. 
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It has been shown recently that intermittency of the Gledzer Ohkitani Yamada (GOY) shell 
model of turbulence has to be related to singular structures whose dynamics in the inertial range 
includes interactions with a background of fluctuations. In this paper we propose a statistical theory 
of these objects by modelling the incoherent background as a Gaussian white-noise forcing of small 
strength T. A general scheme is developed for constructing instantons in spatially discrete dynamical 
systems and the Cramer function governing the probability distribution of effective singularities of 
exponent z is computed up to first order in a semiclassical expansion in powers of V. The resulting 
predictions are compared with the statistics of coherent structures deduced from full simulations of 
the GOY model at very high Reynolds numbers. 



I. INTRODUCTION 

Are structures (sheets or filaments of vorticity) a vital ingredient of intermittency in 3D-incompressible turbulence? 
To date, this important question remains open [1], and an answer starting from first principle, i.e., from a controlled 
approximation to the Navier-Stokes equations, seems over the horizon. The new understanding of the anomalous 
scaling in the Kraichnan's model of passive advection [2], based on the identification of zero modes in the homogeneous 
Hopf equations for equal-time correlators, has in particular strengthened the belief that field-theoretical methods would 
eventually be able to capture the full statistics of turbulent flows without an explicit account of structures. 

Interestingly enough, the relative interplay between coherent ordered structures and incoherent turbulent fluctua- 
tions turns out to be a subtle matter already in the restricted framework of the so-called shell models of turbulence 
[3]. It was noticed very soon [4] that elementary bricks of intermittency in those deterministic ID-cascade models 
could be pulses or bursts of activity growing in an almost self-similar way as they move from large to small scales. 
However, genuine dynamically stable self-similar solutions of the equations of motion in the inertial range display 
an unique scaling exponent (to be denoted below as Zq), provided they are localized in k-space (which, in the shell 
model approach, reduces to a discrete set of wave numbers k„ = 2^ n ~ x \ where the shell index n goes from 1 to oo). 
Furthermore, the exponent zq, giving the logarithmic slope of the velocity gradient spectrum left in the trail of the 
pulse, happens to be rather close to the Kolmogorov value 2/3 (zq = 0.72) in the case of the Gledzer-Okhitani- Yamada 
(GOY) model, in the range of parameters where it reproduces at best the multiscaling properties of real turbulent 
flows. 

In Ref. [5] , the role played by the interaction of pulses with the rest of the flow in producing more singular events 
was unravelled, and a two-fluid picture was introduced, where coherent structures form in and propagate into a 
featureless random background. Our goal in this paper is to elevate this still rather qualitative proposal to the rank 
of a semi-quantitative theory and to test its predictive power about intermittency in the GOY model. We shall 
assume that turbulent fluctuations on the shells downstream the pulse, i.e., small scales, act on the coherent part of 
the flow as a random, white-in-time, Gaussian forcing and ask whether the inviscid stochastic extension of the GOY 
model obtained in this way is able to reproduce the statistics of strong deviations of the full turbulent system in the 
inertial range. There is a priori quite a lot of freedom in the parameterisation of the forcing. Therefore, in order to 
keep things as simple as possible, we bind ourselves to use a single adjustable parameter (hereafter noted T), which 
measures the level of noise. We consider the semiclassical limit T <C 1 of these systems and study the statistics of 
singular structures appearing in this regime. 

Semiclassical (or instanton) techniques are well suited to capture large and rare excursions of fluctuating fields [7] . 
As such, they have gained recently a renewal of interest in the field of turbulence and have already led to noteworthy 
results in the context of Burger's turbulence [9,10], and of the Kraichnan's model of passive scalar advection [11,12]. 
One usually starts from a path integral representation of high order structure functions and uses a saddle point 
approximation to determine the coupled field-force configurations contributing mostly to those quantities. The nature 
of the statistical object to be computed imposes precise boundary conditions on the physical field and the random 
force (respectively at large and small scales, where the cascade processes start and end). Instantons, which in the 
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inertial range often reduce to a self-similar collapse along some spatial dimensions, are eventually selected by a delicate 
matching procedure at the two boundaries. In shell models we are dealing with an intrinsically discrete lattice of 
logarithmic scales. As a consequence, the analytic computation of instantons is completely out of reach in the inertial 
range, not to speak about the matching on both sides of the cascade. To circumvent this difficulty, we shall focus 
on the probability distribution function (pdf) of scaling exponents after n cascade steps, P n (z), and argue that, in 
the semiclassical limit r « 1, this pdf builds up from the neighborhood of a single self-similar instanton (of scaling 
exponent z) that dynamic stability considerations will help us to construct numerically. In order to get non trivial 
physics, it turns out to be necessary to perform the semiclassical expansion of — lim„_ i . 00 ^ In P n (z) (the rate of 
rarefaction of singularities of scaling exponent z in the multifractal picture) up to next to leading order in powers 
of T. This can be achieved via a summation over quadratic fluctuations around the instantons, once the proper set 
of boundary conditions for the corresponding trajectories in configuration space has been defined. We shall show in 
details how to carry out this program and end up with a prediction for P n (z) lending itself to a straight confrontation 
with the pdf of effective scaling exponents of coherents events that can be extracted from simulations of the GOY 
model at very high Reynolds numbers. Although our interest lies primarily in gaining a better understanding of 
intermittency in the framework of shell models of turbulence, the emphasis will be put in this paper on the technical 
aspects of the method that we had to develop for computing instantons. We believe that this method is general 
enough to find applications in other contexts or physical problems, like for instance the motion of complex objects or 
excitations on ID-lattices in the presence of a co-moving random environment. 

The paper is organized as follows. In Section II, we define the stochastic extensions of the GOY model that we 
shall study. In Section III the equations of motion for instantons will be derived using the well-known Martin Siggia 
Rose path integral representation of probability distribution functions for stochastic dynamical systems. Section IV 
is devoted to the computation of self-similar extremal trajectories, with the theoretical considerations underlying the 
solution explained in Subsection IV A and its practical implementing, together with the results, exposed in Subsection 
IV B. The important effect of quadratic fluctuations and the rather heavy formal work behind their computation are 
discussed in Section V. The comparison of the results issuing from the instanton approach with numerical data on 
the statistics of coherent structures in the genuine GOY model is given in Section VI. We conclude in Section VII, in 
particular as to the relevance of a two-fluid description of intermittency in shell models of turbulence. 



II. DEFINITION OF THE STOCHASTIC DYNAMICAL SYSTEM 



Equations of motion for the GOY model in the inertial range read : 

^ = Q 2 (l - + eb* n _,b* n+1 - Q~ 2 b* n+1 b* n+2 , (2.1) 

where the complex variable b n = k n u n should be understood as the Fourier component of the gradient velocity field 
at wavenumber k n = Q n and the integer n runs from to +00. Throughout this paper, usual values of parameters 
e = 0.5 and Q = 2 will be assumed. It is convenient to cast (2.1) in a vectorial form 

§ = N|b], (2.2) 

where the infinite-dimensional vector b is built up from the b n 's, while the rfi 1 component of the nonlinear kernel N[b] 
is given by the right hand side of (2.1). It is worth noting at this point that b* • b = X^^Lo l^l 2 P^ a y s dimensionally 
the role of enstrophy in real flow and that the inverse square root of this quantity sets the order of magnitude of the 
smallest time scale on the shell lattice. 



Since quadratic nonlinearities lead generically to finite time singularities, it is very useful to introduce a desingu- 
larizing time variable r related to the physical time t by the differential law 

§ = (b*-b) 1 / 2 . (2.3) 



This turns (2.2) into 

db _ N[b] 



dr (b*-b)!/2' 



(2.4) 
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where both sides of the equation have the same scaling dimension in the field b, which shows that an infinite "time" 
is now required to form a singularity by travelling across the whole shell axis. 

^,From previous work [13], we know that every initial condition of finite enstrophy, when evolving under dynamics 
(2.4), eventually organizes itself in a soliton-like pulse, moving from large to small scales at a constant speed with 
an exponential growth of its amplitude. The asymptotic state is unique, up to trivial phase symmetries of the GOY 
model [14], time translations or multiplicative rescaling of the field b, which all leave the equation of motion (2.4) 
invariant. We may restrict our attention without loss of generality to the case where the phase pattern along the shell 
axis does not break into a three-sublattice structure. The asymptotic Floquet state, to be noted henceforth b°(r), is 
then purely real and such that : 

b°n+i(T + T ) = exp(A T )b° n (r). (2.5) 

The period To is the "time" needed for the center of the pulse to go from shell n to shell n + 1, while the (positive) 
Lyapunov exponent Aq controls its growth. Both quantities To and Aq, are dynamically selected in an unique way. 
The scaling exponent zo associated to the pulse (fixing in particular the logarithmic slope of the spectrum left in its 
trail) can be extracted from the identity Q z ° = exp(A T ). Its value turns out to be 0.72 in the case of the GOY 
model for the choice of parameters stated before. 

We turn now to the stochastic models, that we are interested in solving by the instanton method. Their physical 
motivation has been explained in Ref . [5] : we assume that pulses parameterize adequately singular (and temporally 
coherent) structures in shell models but that the deterministic dynamics (2.2) should be enlarged towards a stochastic 
one, in order to describe the interaction of a given pulse with incoherent fluctuations produced by the relaxation of 
the trails left by its predecessors. We are therefore led to consider the following extension of the original inviscid GOY 
model : 

^ = N[b] + Vf (b* • b) 3 / 4 B[C]77 , (2.6) 
where 7] is a Gaussian noise, delta-correlated in time and shell index, whose correlations read : 

(V*n(t)Vn>(t'))=5 nn ,5(t-t'). (2.7) 

The various factors coming in front of rj in (2.6) have the following meaning : the number T fixes the relative strength 
of incoherent fluctuations with respect to coherent ones and we shall be interested in the semiclassical limit of small 
T amenable to semi-analytic treatment. As will be clearer in a while, the overall scale factor (b* • b) 3 / 4 is there to 
keep noise relevant all along the cascade, thereby preserving scale invariance. Finally the matrix -B[C], of zero scaling 
dimension in the field b since it depends only on the unit vector C = , b , may be used either to introduce spatial 

v( b * b ) 

correlations of noise (along the shell axis) or to localize its action with respect to the instantaneous position of the 
pulse. Although the formalism to be developed in this paper can deal with the most general situation, we restricted 
ourselves in practical investigations to diagonal matrices -B[C], just playing with the degree of localization of noise. 
Results will be presented for three rather emblematic choices of B : (i) B nn = 1, which describes a completely 
delocalized noise; (ii) B nn = C*_ 2 C*_ 1 which keeps some flavour of the original GOY dynamics and makes noise 
active just at the leading edge of the pulse; and finally (iii) B nn = \C n -5\ 2 + \C n -4\ 2 , which removes the action of 
noise further away from the center of the pulse. We must emphasize that these particular choices were not dictated 
by rigorous considerations on the underlying dynamics of the GOY model, but rather used to scan the variety of 
behaviours which may be expected from such stochastic dynamical systems. It should be noted that the strucutre of 
the matrix B is not constrained by any conservation law, since the coherent part of the flow does not form anymore 
a closed system, even in the inertial range, in our two-fluid description. Finally, to simplify the following analysis, 
we are going to restrict the fields b and r) in (2.6) and (2.7) to being real-valued vectors and neglect the effect of 
imaginary fluctuations. This is certainly not a serious restriction as for the instantons themselves, which are expected 
to be, like the self-similar deterministic solution described above, purely real, up to trivial phase symmetries of the 
GOY model. It can also be remarked that the model (ii) (which will be found later on to give the more convincing 
results) does not require a complex noise, since the phase degrees of freedom have already been incorporated in the 
definition of the matrix B in that case. 

While the deterministic dynamics (2.2) selects a single self-similar solution blowing up in finite time with scaling 
exponent zq, the presence of noise in (2.6) allows for a continuum of scaling exponents, even in the manifold of 
normalizable fields b. In the small noise (or semiclassical) limit r C 1, the probability density of developing an 
effective growth exponent z after n ^> 1 cascade steps will take the form : 
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P n (z)~Vn exp 

where sq(z) is the action per unit cascade step of the self-similar extremalsolution of scaling exponent z of optimal 
bare Gaussian weight (or instanton), and Si(z) measures, to lowest order in T, how the basin of attraction of the 
instanton in phase space evolves with the number of cascade steps. Note that the argument of the exponential in 
(2.8), — ( SO p^ +Si(z)), is nothing but the Cramer function introduced in the theory of large deviations, which governs 
the rate of rarefaction of singularities in the multifractal picture [15]. We will show in this paper how to compute 
in a clean way both quantities sq(z) and s\(z). Before doing this, we must carefully handle problems related to the 
time discretization of the stochastic equation (2.6) since a consistent treatment of them is necessary to get the right 
expression of the first correction Si(z). We shall adopt the view that the initial stochastic equation (2.6) is to be 
understood in the Stratonovich sense [16]. However, in the path integral formulation of stochastic dynamical systems 
that we shall heavily use in the following, it is much simpler to work with the Ito prescription which, in the limit of 
small time steps, amounts integrating (2.6) within a basic Euler scheme with all b-dependent quantities in the r.h.s. 
estimated at the prepoint. When switching to the Ito discretization recipe, the stochastic equation has to be changed 
into : 

^ = N r [b] + Vf (b ■ hfl A B[C] n , (2.9) 

where the new kernel Nr[b] differs from N[b] by the addition of the so-called Ito drift-term. We give, for the sake of 
completeness, the resulting expression of the n th component of Nr : 

N rn [h] = N n [b] + A[( b • b)^ 4 B nj ](b ■ h)V 4 B kj . (2.10) 

At this point, we may write down the discrete analog of (2.3) as Ar« = r, + i — t, = (b, - bj) 1 / 2 (ii+i —U) (where i is the 
time index) and redefine the noise as r] i — > ^ = V / T(b, • b i )~ 1 / i r} i . This leads to the following stochastic extension 
of (2.4) which will be the starting point of our formal analysis : 

with 

(Ut)Zu' (t')> = TS nn ,S( T - r') . (2.12) 



sq{z) 



+ si{z) 



(2.8) 



III. EXTREMAL TRAJECTORIES FROM PATH INTEGRAL FORMULATION 

Statistics of classical fields in the presence of random forces can be examined with the help of field theoretical 
techniques formulated in [17]. In particular, the probability to go from point b, n at time r = to point b/ at time 
Tf may be written as a path integral : 

P(b in ,0;bf,T f ) = JvbVpexp-S[b,p], (3.1) 

where 5[b,p] is an effective action to be defined below, p an auxiliary field conjugated to the physical one b and 
VbVp stands for : 

11 O^d ■ ^- Z > 



W 11 (27T) 



In the last equation, the time interval r/ was divided into TV subintervals of length At = ^ (with b in = b and 
bf = bjv) and the number of shells was set to a finite value d, in order to give a clear meaning to the measure. For 
the problem of interest (2.11), the effective action S takes the form : 

i =N ~ 1 / TV fK 1 \ r 

5[b,p]= J2 ( b ^+i ~ b ^ ~ At (b . .b.)i/2 j + 2 {hi ' hi) Pi ' B l h iV B l h i]Pi ■ ( 3 - 3 ) 
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or, in the continuum limit : 

S[b, p] = £ drip. - (b Nr b [ ) b 1 ] /2 ) + |(b • b) p • B[b] *B[b]p . (3.4) 

The last term in (3.4), quadratic in p, appears as a result of averaging over the Gaussian noise £, while the first one, 
linear in p, would still be there in the absence of noise as a formal way of enforcing the deterministic equation of 
motion of b. S[b, p] will be referred to in the following as the Martin-Siggia-Rose (MSR) action. 

Rescaling the auxiliary field p as p'/T puts an overall large factor 1/T in front of the effective action and opens 
the way to a saddle point approximation to the path integral (3.1). Extremization of the action with respect to the 
configurations of both fields b and p between times and r/, for fixed endpoints, leads in a straightforward way to 
the following set of coupled equations defining extremal trajectories : 

rW 1 

— = - t M0--d b [(b-b)0-B t B0]. (3.6) 

(IT 2 

In the above equation, we set p' = — i 6 and M. is the Jacobian matrix of the kernel ^^|^j 2 : M. = 

- ™J£- As usual, equations (3.5) and (3.6) inherit a canonical structure 
(b • b) L ' z (b • b) A l z 



db _ &H 

dr~ ~ 1)6 ' 
dO _ dH 
~dr~ ~ ~~db 



(3.7) 
(3.8) 



where the Hamiltonian T-L reads 



n= ^w^ + l {h - h){e - Btm) - (3 - 9) 

Since H is not explicitly time-dependent, we conclude that its value, to be called the pseudo-energy in the sequel, is 
conserved along any extremal trajectory. The action S[b,0] may be rewritten in terms of % as 



S[b,0] 



r *(•■£-")■ <3io) 



from which it is seen that the further requirement that the trajectory be extremal with respect to time reparametriza- 
tion leads to the condition of vanishing pseudo-energy H = 0. Noting that each term of H in (3.9) has the same 
scaling dimension in b and 6 (either 1 or 2), one finds that 

which shows that the overlap b -0 between the physical and auxiliary fields is also conserved, together with 7i, along an 
extremal trajectory. This property reflects the scaling invariance of the stochastic cascade processes we have in mind. 
We should at this point insist on the fact that, in contrast to instantons in the framework of equilibrium statistical 
mechanics or quantum mechanics, equations for extremal trajectories in stochastic dynamical systems describe the 
real motion of the physical field in a particular "optimal" realization of the noise. The comparison of equations (2.11) 
and (3.5) shows indeed that the following relation holds between £ and 9 : 

$, = {bb) 1 / 2t Be. (3.12) 

Like their deterministic parent (2.4), the equations of motion (3.5) and (3.6) sustain formally traveling wave-like 
solutions, such that 

b n+1 (T + T) = expATb n ( T ), (3.13) 

9 n+1 (T + T) = exp-AT9 n (T), (3.14) 
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whose scaling exponent z = is expected to be now related to the overlap /Ji = b • 6 (with z = Zq for /Ji = 0, in the 
absence of noise). However there is little hope to find these solutions by a direct forward in time integration of (3.5) 
and (3.6), as could be done successfully for equation (2.4). This is because the auxiliary field 6 intrinsically propagates 
"backward" in time, as is clear from the discretized version of (3.6) (deduced from the extremization of (3.3)). In the 
present problem, we have observed numerically that regular Floquet states emerge as dynamical attractors of (3.5) 
and (3.6) only for rather high values of /ii (otherwise the system evolves in a chaotic manner). They form a branch 
of solutions definitely distinct from the one to be obtained in the next Section and correspond presumably to local 
maxima of the action rather than the local minima of interest to us. 



IV. AN ITERATIVE METHOD FOR COMPUTING SELF-SIMILAR INSTANTONS 



A. Theory 



The previous considerations suggest that equations (3.5) and (3.6) should not be treated on the same footing. 
The careful examination of physical properties that instantons should possess will give us keys for computing them. 
Assume for a while that a solution has been found, obeying to (3.13) and (3.14). We note b°(r) and £°(t) the 
corresponding configurations of b and £. The linearization of the equation of motion (3.5) at fixed noise leads to the 
following evolution of fluctuations Sb of b around b° : 

jUb = LSb=MSb + (Sb • 9b) ((b • b) 1 / 2 ^ . (4.1) 

The periodicity properties of the linear operator L ensure that the fluctuations of b may be decomposed on a complete 
set of eigendirections ^^(r) evolving according to (4.1) and such that 

* ir (r + T) =e^ T T+i* if .(r), (4.2) 

where 7+i denotes translation by one unit in the right direction along the shell lattice. In practice we shall have to 
work with a finite number of shells d and, in order to get rid of boundary effects, it will be necessary to fully periodize 
the shell lattice : the index i then runs between 1 and d and the translation operator is easy to represent as a matrix. 
Formally, the ^i r 's can be determined at time r = by diagonalizing the Floquet operator : 

U T = T-iexp f Cdr, (4.3) 
Jo 

where exp is a chronologically time ordered product (initial time on the right). One observes that b° satisfies (4.2) 
with a time averaged Lyapunov exponent a = A. 



We claim now that every initial condition b(0) evolving in the configuration of noise £ should be attracted towards 
the instantonic trajectory. If it were not true, some perturbations would be able to grow in the comoving frame of the 
pulse, thereby generating scaling exponents larger than z at no cost in the action, in contradiction with the hypothesis 
that the optimal realization of a singularity of exponent z has been found. This strong criterion of dynamic stability 
is another way of stating that the Cramer function should be insensitive to the details of the production of pulses in 
the forcing range. It implies that A is an upper bound for the real part of the <r,'s. Arranging the eigendirections 
in order of decreasing Re<Ti, we are therefore led to identify b°(r) with H>i r (r). In the case of zero noise where we 
recover the deterministic solution of Section II (with £° = in both (3.5) and (4.1)), the time derivative is also 
solution of (4.1) with the same Lyapunov exponent as b°, a = A (= A in this case). In this limit we would naturally 
define ^2t(t) as This property is lost in the more general situation of a non vanishing noise, because £° is not 

time- invariant. What remains true however is the fact that b° and still span the set of "coherent" fluctuations 
which do not affect the shape of the pulse but modify its height and position. 

By turning now our attention to the linear dynamics dual to (4.1) we shall come close to the equation (3.6). Let 
us indeed consider the equation of motion 

£ = (4-4) 
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where in order to limit the proliferation of symbols, we keep the same notation 6 for the new auxiliary field, although 
it is only in particular circumstances, to be clarified below, related to the of equations (3.5) and (3.6). From (4.1) 
we get : 

M=-t Md _ dh £ (b . b) l/2 g | _ (4 5) 

The dual dynamics enables one to construct a basis of left eigenvectors *f? u (r) (with 1 < i < d) satisfying 

*u(t + T) =e-^ T T + i* i /(r), (4.6) 
as well as the following orthogonality conditions with the members of the first basis : 

*«(T)-* jr (T)=<f«, (4.7) 

at every time. The vectors ^a(0) are determined by diagonalizing the adjoint Floquet operator 

i c/ T = exp f '£drT+i, (4.8) 
Jo 

and enforcing the normalization condition (4.6) at time r = (exp is now an anti-chronologically time ordered 
product). We may note at this point that the first left eigenvector *f?u is, in the generic case of non-zero noise, 
the only one to display the scaling behaviour anticipated for 0° according to (3.14), since its Lyapunov exponent 
equals — o\ = —A. We conclude that the auxiliary equation (3.6) in the restricted manifold of self-similar solutions is 
tantamount to the relation : 

0°(t)=pi*i«(t), (4.9) 

where the multiplicative constant fii is nothing but the overlap 0° ■ b° (= ^/i*u(r) • *i r (r) = Hi), which was shown 
before to be indeed a conserved quantity. This claim is further confirmed by rewriting the original equation (3.6) as 

rlf) 1 

— = - t M0-(6-B t B6)h-(bh)6.-[(d b B) t B + B(d b t B)}6 

UT 2 

= - l M 6-(6-B 'B0)b - (b • b) 6>.(<9 b B) *B6> . 
Putting back superscripts and reintroducing £° by using (3.12), we arrive at 

./flO i o 

d JL = _t M e o _ ( 0O . B ^ ) _b_ _ (b o . b o )1/2 . B€ oj (4 

which shows that 0° obeys the dual dynamics defined by equation (4.4). 

Having interpreted (3.6) as a condition of self-consistency for the conjugate momentum 6 expressed by (4.9), we 
could contemplate the following Newton- like procedure for catching numerically self-similar instantons. First make 
a guess for £ in the form of a traveling wave (£ m (r + T) = 7+i £ m ( T ))i integrate (3.6) forward in time in order to 
determine the asymptotic traveling state reached by b in the prescribed configuration of the noise. Then compute 
S&i/ from the diagonalization of 1 Ut (or from running (4.4) backward in time in order to let emerge the eigendirection 
of lowest growth rate), employ (4.9) for producing a new configuration of (and thereby £) and iterate this loop 
many times at a fixed value of the overlap /ii until convergence is achieved. However two major difficulties call for an 
improvement of the method : they both have to do with the stability of the trajectory upon time reparametrization. 
First we do not know the speed (or the inverse period T _1 ) of the final traveling wave which must carry together b and 
£. Therefore, when performing the first step of the iterative loop, we must allow continuous time reparametrization of 
our Ansatz for the noise in order to fine tune the speeds of the two pulses formed by b and £ and let both terms in the 
right hand side of (3.5) be always relevant. It will be explained in the next Subsection how this goal can be achieved in 
practice. The second difficulty is much more serious than the preceding one and in the way of getting around it resides 
perhaps the most tricky part of this work. The point is that traveling wave solutions to (3.5) and (3.6) may perfectly 
have a non zero pseudo-energy 7i, while we are looking for the particular ones with V. = 0. We shall be able to fulfill 
asymptotically the two conditions h.O = Hi and H = 0, if and only if our iterative guess for 6 is constructed within a 
two-dimensional space rather than a unidimensional one as in the naive proposal made above. For this purpose, we 



7 



are going to embed the linearized dynamics (4.1) into a new one which admits the time translation mode as a true 
eigenstate, of the same growth factor A as b°, restoring thereby the symmetries present in the absence of noise. We 
shall do that in the most economical way, from both formal and numerical points of view, by substituting <J> 2 r = 
to *2r, i-e., the eigendirection along which the fluctuations of b around b° are the less stable. 

The left eigenvector *&2i is first rescaled as $2; = ir.^^L — (which makes sense as long as \&2!-3»2r ^ 0, a condition 

always found to be satisfied in practice), so that &21 has a unit overlap with <& 2r , while being orthogonal to all other 
right eigenvectors \&j r with i ^ 2 (which will be noted from now on). One then considers the modified linear 
dynamics : 

^ = CSb = C6b + (b° • b°) 1/2 (* 2/ • 6b) B ^y- . (4.11) 
dr dr 

It is easily checked that <£ 2r = ^7- obeys (4.11), since equation (3.5) yields upon time derivation : 

^=£* 2r + ( b°.b ) 1 / 2 i?^=£* 2f ,. (4.12) 
dr dr 

It is also trivially seen that the other vectors = Slv for i ^ 2 keep the same evolution under (4.11) as under (4.1). 
The dual dynamics reads now : 

^ = - l t0 = - f C9 - (b° • b ) 1 / 2 (0 ■ B ^)* 2/ . (4.13) 
dr dr 

It leads to a new family of left eigenvectors <!>,;, dual to the direct basis, whose second member has been defined 
above and the others relate to their original counterparts as 

*iZ = *iZ " • *2r)*2/ • (4.14) 

Although the &u's were introduced as a rather formal trick, it should be emphasized that $1; and $ 2 ; have an 
appealing physical meaning. Parameterizing a perturbed trajectory for b as b = e S]nb ( T ) b°(r + St(t)) + Sbi nc , where 
the "incoherent" part of fluctuations Sb inc is bound to be a linear superposition of the less dangerous modes <£j r for 
i > 3, one has indeed, to linear order in Sb, 

5\nb = #u ■ Sb , (4.15) 
ST = $ 2 r$b. (4.16) 

These two relations will be useful in the computation of quadratic fluctuations to be presented in Section V. They 
show that by projecting out the multidimensional fluctuation field Sb onto the two vectors $1; and $ 2 ;, one has 
access to the most relevant part of it affecting respectively the amplitude and the time delay of the pulse constituting 
the instanton. In terms of <l>i/ and <& 2 /, the self-consistency condition (4.9) for 0° together with the requirement of 
zero pseudo-energy H take the following form : 

0°(r) = A*i*nto + Mt)* 2; (t) , (4.17) 

where 

H2(r) = \e° -B^BO = \i a -e ■ (4.18) 

Since from (4.17), /i 2 (r) = * 2r .0°, and from the equation of motion (3.6), $ 2r • 6° = U + • £°, the relation (4.18) 
is just a way of restating % = 0. That (4.9) implies (4.17) results from the general link between \l>i/ and *i/ (see 
(4.14)). The reverse is true only under the supplementary condition of constant % or /< 2 (r) = C te + • £°, which is 
guaranteed by (4.18). It is proven by checking that in that case 0°(t), as given by (4.17), obeys, as it should, (4.4) : 

dr dr 

= -tjCO - (b° • b°) 1/2 (0° • B ^)* 2; + ^ * 2; 

dr dr 

dr 2 dr 



8 



The great advantage of (4.17) and (4.18) with respect to (4.9) is that this couple of equations lends itself to iterative 
procedures leading inexorably to a fixed point of zero pseudo-energy, a task seemingly out of reach before. There is 
some unavoidable arbitrariness in the construction proposed here, concerning in particular the definition of the vector 
about which the reader may feel a little uncomfortable. We suspect that these unwanted features do not affect 
the final results since the original equations to be solved as well as the corresponding conserved quantities all have a 
clear mathematical definition where <l>i/ and <1>2/ merge together into 



B. Practical implementing and results 

The action density Sq(z) could be computed successfully for the three stochastic models defined in Section II using 
the iterative scheme outlined before. The shell lattice was first mapped onto a circle of d sites, with d typically ranging 
between 20 and 30. Finite size effects turn out to be completely negligible at such lengths of the chain, due to the 
strongly localized structure of the instantons. To start the computation, we make a guess for both the unit field 
^( T ) = (b b)V^(r) anc ^ ^ e n °i se £( T ) cane d henceforth C m (r) and £ m (r). They are such that 

& n + + T<n) = c < n(T)j ^„ +i(r + T<n) = £, (t)j (4 1Q) 

and 

CT +d {r) = Ci n (r), C +d (r) = C(r). (4.20) 

Furthermore, the noise is normalized in such a way that the overlap b • 6 = C • B _1 [C]^, takes on a prescribed value 
Hi held as a control parameter during all the steps of the computation. A possible and convenient choice would 
be for instance C m (r) = C°(r), where C°(r) is the deterministic solution of scaling exponent zo, and £*™(r) = 
fii(b° ■ b°) 1 / 2 B[C°]*? / (r), where *i/(t) is the left eigenvector dual to b°(r), i.e., in the absence of noise. Fig. 1 
shows how both vectors and like at a given time. Their shapes will in fact little evolve as we let the 

scaling exponent z depart from Zq. 

In order to allow time reparametrization of the trajectory, we first get an estimate of the instantaneous position of 
the pulse along the shell axis in our trial configuration by computing the following quantity 

n»(r) = f «M(C(r)) 2 . (4.21) 

n=0 

The notation n[d] recalls that, due to cyclic boundary conditions, the shell index n is now only defined modulo d and 
that in practice a continuous determination of this integer should be adopted close to the center of the pulse which 
contributes mostly to the right hand side of (4.21). One has by construction n m (r + T m ) = n m (r) + 1 (mod d). 
Having recorded n m (r) and £ m (r) during a whole period T m , we integrate forward in time the nonlinear evolution 
equation for C deduced from (2.11) by projecting out the longitudinal part of its right hand side : 

^ = {N[C](r) +B[C](r)e n (r')} ± , (4.22) 

Note that the subindex T disappeared from the nonlinear kernel Nr because the Ito-drift term being linear in T does 
not matter in the computation of the action to leading order (we shall see in the next Section how to handle it to next 
to leading order). The most salient feature of (4.22) is that the noise configuration is evaluated in relation not to the 
time r but rather to the actual instantaneous position of the pulse. This means that the time t'(t) is automatically 
delayed or advanced with respect to r, according to the recipe 

n in (T')=n(T) or r' = (b'T'Nt)]. (4.23) 

After integrating (4.22) long enough, a new traveling wave state C oui (T), £°"'(t) = ^ m (r'(r)) will usually emerge, of 
period T out possibly different from T m and averaged growth factor 

A out = ^ J + (N[C oui ] + B[C out ]£ out ) .C out dr. (4.24) 

The vectors *i/(r), *2/(t) are then identified as the two eigenvectors of *?7t°»* (defined in (4.8)) of smallest (real 
negative) Lyapunov coefficient and the corresponding *i/(r), *2/(t) constructed as linear combinations of them 
obeying for all times the following relations 
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Su-b "' =*«.—— = 1 
ar 

* 2/ -b™' =*!,.__ =0 
ar 

Finally, the trial noise configuration is renewed as £ in (r) = (b oui • b™') 1 / 2 '5[C°"']0 in with 

irl (r) = /ii*i/(r) + /i 2 (r)* 2/ (r), (4.25) 
where /i 2 (r) is determined upon imposing the condition of zero pseudo-energy on the trial solution (0 m (r), b otl '(r)) 

gin . N [ C o«i] + 1 ^out ^out^gin . g^qout^ t^out^gin = q ^ 26 ^ 

After setting b m (r) = b°"'(r), we are ready to repeat the operations described above as many times as needed 
until a fixed point of the transformation (such that b oui (r) = b m (r) and £ m (r) = £ ou< (t)) is reached, which solves 
the problem. The good stability properties of the algorithm, as well as its iterative character, authorize a rather 
unsophisticated handling of issues raised by the time discretization. As required by the Ito- convention, the equation 
of motion (4.22) was integrated using a first-order Euler scheme with a time step Ar = about 350 times smaller 
than the period. The time r' was approximated as the multiple of Ar making the relation (4.23) best satisfied. 
Similarly no higher order interpolation scheme was devised for estimating with accuracy the output period T° ut : 
it was again simply approximated as the multiple of Ar making periodicity conditions (4.19) best satisfied for the 
output. However the time step Ar was changed at each iteration of the loop so to maintain the time resolution N 
constant. The efficiency of the method was greatly improved, when seeking solutions of exponents z far from zq, by 
increasing z gradually (through an increase of the control parameter /ii ) using as first guess solutions of lower but 
close scaling exponent z'. In this way convergence toward satisfactory self-similar solutions (of exponent z varying by 
less than 10~ 5 under iteration and pseudo-energy H < 10~ 5 ) was attained in no more than 20 iterations. 



We turn now to the presentation of our results. Fig. 2 shows the action density sq(z) for the three models (i), (ii) 
and (iii) defined in Section II. Values of T were adjusted in order to provide the same curvature of sq(z)/T at the 
bottom of the curves, reached evidently at z = zq. We see that the variations of the zeroth-order action get sharper 
on the z > zq side (the only one displayed in Fig. 2), as one goes from model (i) to model (iii). Figures 3, 4 and 5, 
referring respectively to models (i), (ii) and (iii), show the normalized coherent field C and the random force B£ at 
increasing values of z (0.75, 0.85 and 0.95). In all cases the random force is found to be negative at the leading edge 
of the pulse, in agreement with the physical picture advocated in [5] : growth can be enhanced only by frustrating 
the energy transfer processes. For model (iii), the coherent field itself gets negative at the forefront : noise in that 
case just helps to prepare the system in an initial condition consisting of a pulse and a negative well in front of it, 
which then collide. An interesting upshot of our computations is that models like (ii) or (iii) involving only a local 
coupling between b and £ escape the disaster met in the framework of model (i), namely a cross-over toward an 
asymptotic linear growth of Sq(z) with z, already perceptible in Fig. 2. Such a behaviour forbids the existence of 
velocity moments at arbitrary orders and is, thus, clearly undesirable. It turns out that the whole shape of sq(z) for 
model (i) can be pretty well understood from an adiabatic approximation, which is carried out in the Appendix A, 
where solutions of arbitrary scaling exponents are constructed using adequate time reparametrizations and dilations 
of the deterministic solution of scaling exponent zq ■ The validity of this approximation for the model (i) is somehow 
obvious from Fig. 3, where it can be checked that instantons keep indeed almost the same shape, even for quite 
sizable variations of z. Its failure in models (ii) and (iii) is probably due to too strong deformations of the solutions 
as z increases, again suggested by Figs. 4 and 5. The full non-linear treatment of the problem proposed in this paper 
was however necessary to reach this quite fortunate conclusion. 



V. THE EFFECT OF QUADRATIC FLUCTUATIONS 



A. Formal considerations 



In order to compute the first order (in Y) correction to the action per unit cascade step (si(z) in the expression 
(2.8)) of the density of probability P n (z), we have to expand the MSR action up to quadratic order in fluctuations Sb 
around the extremal trajectory b°(r) of scaling exponent z and then sum over them in a way which will be explained 
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below. Since typical fluctuating paths are not differentiable but rather behave as Wiener paths with derivatives of 
the order i, we shall stick to time-discretized expressions in all the following manipulations of the path integral. For 
the sake of clarity, the superscript referring to the extremal trajectory in the previous Section will be taken away, 
whereas bj, 6i will be short-hand notations for b°(rj = i At), 0°(Tj = i At), where St is the (small) time step used 
in the discretization. 

We start from (3.1) and the representation (3.3) of the effective action 5[b,p], expand it to quadratic order in both 
fluctuations Sp and Sb and then integrate out the fluctuations of the auxiliary field. To begin, the time t/ during 
which we let the system evolve will be equal to the time nT needed by the ideal instanton to perform n steps along the 
shell axis. The ideal initial and final configurations b in = b and b/ = bjy describe then a pulse centered successively 
around the shells of index and n. To quadratic order in deviations from the extremal trajectory, the probability of 
joining the perturbed endpoints bj n = bo + Sbo and b/ = bjv + 5bjvin the time t/ take the following expression : 



P(b + Sb , 0; b N + Sb N , Tf ) = e - n '°W T j VSb exp (- J5[*b] - S 2 S[Sb}) , 
where the measure of integration is defined as 



VSb 



1 



1 



1 



2tiTAt J (b • b ) d / 2 | detS 



JV-l 

n 



ddbj 



(bi-bi)^ |detBi 



the linear variation SS reduces to the boundary term 

1 



SS[db] = -{On-! ■ Sb N - ■ db } . 



and the quadratic one reads 



<5 2 S[<5b] 



2T 



JV-l 

E 

8=0 



AiSbi) 



+ Sbi.ViSbi 



The drift- and potential- terms in eq.(5.4) are found to be given by the following relations 

AiSbi = MiSbi + (Sbi ■ db,) ((hi ■ bi)Bi %0i) , 

and 

1 



i .(8b i -d bi )M i 8b i --0 i .{8b i -d bi ) 2 ((bi-b^Bi'BiOi) , 



(5.1) 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



(5.6) 



Sbi.Vidbt 

where M is, as in previous Sections, the Jacobian matrix of the nonlinear kernel ^^72 • 

Our task is to perform explicitly the integration over fluctuations Sbi, • • • , <$bjv_i at intermediate steps in the path 
integral (5.1). First it will be convenient to get rid of the anisotropy of the "mass" tensor acting in the kinetic term 
of S 2 S [21] by switching to normalized fluctuating fields defined as 



Sbi = (bi ■ bif/ 2 Bihi 



In this way the measure in the integral transforms into 

Vb. 1 



(bo-bo)^ 2 |detfl„| 

i=N-l 



T>h with 



2ttFAt 



n dhi 



(5.7) 



(5. 



In performing this change of variables in S 2 S, we must pay attention to the fact that Sbi + i — Sbi, as well as h, + i — h,, 
are potentially of order At 1 / 2 . One finds that, up to 0(At 3 / 2 ) corrections (negligible in the continuum limit), S 2 S 
becomes : 



S 2 S[h] 



2T 



N-l 

E 

i=0 



1/2 /hj+i 



At 



hi hj + i + hj 

"I" ^ 



A'ihi) 



+ hi.V'hi. 



(5.9) 
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where 



Qi 



1 + 



b»+i ■ bj+i 
b, ■ b 7 



1/2 



/ 77 tjD — 1 



1 + 



b 7 • b, 



1/2 



i+l 



(5.10) 



Vl = (bfb^BiViBi, 



and 



Ci = ~h [ (bi+1 ' b ^+ 1 ) 1/2 ^+ 1 + ( b * ' b <) 1/2 £<] 1 [(b i+ i.b i+1 ) 1 /2B i+1 - (hi ■ bif^Bi 

A few simplifications are now in order. On one hand we expand Qi as Qi = 1 + SQi + 0(At 2 ), where 

/ bj+i ■b i+1 \ 1/2 '.Br 1 + a- 1 £ <+1 _ 
V I) l> y 2 ' i 



(5.11) 



(5.12) 



(5.13) 



(5.14) 



is of order At. Tracking O(Ar) terms in <5 2 S, we may replace Q, by 1 everywhere in (5.9), except for the kinetic 
term _L (h«+i-h.).Q.(h, + i-h,) ^ w j 1 j c j 1 accorc Jing to standard computation rules [18] in path integrals can be reduced to 

W ^ h ' + Ar h '^ — ^ Tr^Qj. On the other hand, exact ways of tracing out quadratic forms like S 2 S are available only once 
they are written in terms of the mid-point field h ' + ^ +h * rather than h, (which amounts going back at this stage to 
the Stratonovich prescription). Again to order At, the substitution of h ' +1 2 +h * to hj in (5.9) is harmless except in the 



cross-product term — f (h i+ i — hj) • -4^h, which gets into — f (h i+ i — h,) • A\ 



, h i+ i +hi At 



Putting things together, we arrive at the following expression for the transition probability in the neighborhood of 
the instanton : 



P(6b in -> db f ,T f ) =| detjjjjfo) | e -s (r t )r e -x 1 (r t )^ { e II - 1 .S bt -e .S bi „} z ^ in ^ h/jT/]j 
where we defined 

1 N-l Ar «-1 

8=0 8=0 

and the reduced path integral Z[hi n — > h/,r/] as 

Z[h,„ -> h/, T f ] = [ " f Vh exp -S 2 S[h] . 



with 



N-l r 



= £ E 



i=0 L 



h i+1 + hj + i + hj 



and 



One easily shows that in the At — > limit 

TrSQi > Tr In 



(5.15) 



(5.16) 



(5.17) 



(5.18) 



(5.19) 



(b A r • b N yl 2 B, 



(b ■ b ) 1 /2 J B 



dATt 
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since the periodicity of the instanton implies that detBjy = detBo, while 



However the ensuing expression for 2i(r/) is not yet complete to order 0(T°). This is because we discarded the 
O(r) Ito drift-term in our computation of extremal trajectories. The small deviations induced by this term may be 
neglected in Ii(t/) and — lnZ[h, n — > h/,r/] which are already first order corrections in a T-expansion but they must 

s (z) 

be taken care of in the zeroth order term. Rather than n — — — , it should read Sr[b r ], where the index T denotes 

a quantity or field evaluated in the presence of the Ito-term specified in (2.10). To first order in T, we can take 
advantage of the extremum property of b r and write down 

Sr[b° r ] - S[h°] = S r [K] - S r [b°] + S T [b°] - S[b°] 

N(b°]-N r [b°] 



c r h 0l cruoi _ 1 f Tf ff o N(b»]-N r [b»] 
Sr[b]-S[b]--y o 6. (b0 . b0)1/2 dr, 



where the last identity just comes from the expression (3.4) of the action. Rearranging things under the assumption 
(satisfied by each of the particular models that we considered) that the entry Bjk of the matrix B involves homogeneous 
monomials of degree I built up only from components C m with m ^ j (one has I = for model (i) and I = 2 for 
models (ii) and (iii) defined in Section II), one finds that Zi(t/) in (5.15) should finally be understood as : 



X 1 (r f ) = ^Ar f + ^^r^.B^ d r. 
2 4 J 



(5.20) 



This preliminary work being done, the discussion will now concentrate on the reduced path integral Z[hi n — > h/, t/]. 
We restrict first our attention to the case of fixed endpoints h in = h/ = 0. The instantons found in the previous 
Sections are physically satisfactory only if the quadratic functional <5 2 S[h] is positive for all {hj} such that h = h/v = 
(we shall see later on that this is not a sufficient condition in the present problem). According to standard results of 
functional analysis [19], positiveness of 5 2 S'[h] is tantamount to the absence of points conjugate to the origin during 
the whole time interval [0,r/]. Recall that the definition of conjugate points goes as follows : d being the dimension 

of the space (here the number of shells) we construct d initial conditions (\i a \h^) such that 

>»$=0. h % ] = At 5 a p + 0(At 2 ) , 

where j3 is a shell index running as a between 1 and d and let them evolve under the Euler-Lagrange equations derived 
from 5 2 S'[h]. The time r, = iAr is said to be conjugate to the origin if the system formed by the d vectors gets 
degenerate there. One can build a matrix Ui such that U"^ = /if Q , in terms of which the initial conditions read 

U o =0, Ui = At + 0(At 2 ), (5.21) 

while Ui+i (for 1 < i < N — 1) is obtained from £7«_i and Ui through a matrix Euler-Lagrange equation, derived from 
(5.18) and conveniently cast into the following form : 

+ - = \v,SLt^ + \V U "tijtB, (5.22) 



where 



Pi = u i+1 -Ui _ Bi u i + u i+1 ^ 



can be seen as a matrix momentum. In the absence of conjugate points, det [/, never vanishes except at the origin 
and one gets the following simple expression for the reduced path integral, provided the At — > limit is ultimately 
taken : 

*„ = „ 1 , = ,,, ] = (i f ) i '*- s i=. (5.24) 
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Details on this result, which may be found in many textbooks on path integrals [20,21], are provided in the Appendix 
B. A very nice feature of the proof of the connection of the positiveness of <5 2 S[h] with the absence of conjugate 
points is that it provides also an efficient way for computing the reduced path integral with h/ arbitrary (but still 
small naturally). Indeed the main idea consists in adding to <S 2 S[h] a boundary term of the form 

1 N1 

S 2 S' = -— ^2 ( h *+i • W i+1 h i+1 - hi ■ Wihi) , 

i=0 

where W is a symmetric matrix at all times and then selecting the right W such that S 2 S + S 2 S' becomes a perfect 
square. In order to achieve this task, W, must be a solution of a matrix Ricatti equation (see again the Appendix 
B for details) which, upon the substitution of a new unknown matrix Ui defined implicitely by the relation (for 
< i < N - 1) : 

^ /TTT J J . TIT T J N Ui + \—Ui Ui+Ui + l 

-(W i U i + W i+1 U i+1 ) = Bi , (5.25) 

Z At Z 

is found to be nothing but the matrix Euler equation (5.22). It follows that the result (5.24) may be extended to the 
case of an arbitrary final configuration h/ as 

Z[h m = - h„r,] = (-L)"" - 7 =L= exp-i.h; ■ W(r f)hf , (5.26) 

where U{tj) and W{tj) are to be computed by letting both matrices evolve from r = to r = r/ according to (5.22) 
and (5.25). There are some subtleties about the choice of initial conditions and proper counting of the number of 
unknowns whose discussion we prefer to relegate in the Appendix B. 



B. A physical definition of s\(z) 

We are now in a good position to compute the next to leading order term s\(z) in the expansion of 
— lim n _s. +00 ^\ogP n (z) in powers of Y. We shall obtain an estimate for P n {z) by summing over all the trajecto- 
ries which lead to the same growth of the pulse as the ideal instanton b°(r) after n cascade steps. In the small V 
limit, all statistically relevant trajectories remain close to b°(r) and we may define unambiguously their "arrival" 
time at the shell of index n as the first time r n such that 

b(r„) =b°(r° =nT) + Sh', (5.27) 

where Sh' reduces to a linear combination of stable "irrelevant" modes $, r (r°) for i > 3. Up to multiplicative factors 
growing at most algebraically with n, we can then write P n {z) as the following integral over the arrival time and the 
position of the endpoint 

P n (z) » J P(b°(0) ^b (T°) + 5b',T n )^(*l/(T°).5bV(* 2 /(T°).5b')dT„^b'. (5.28) 

The density of probability in the right hand side of this expression is known from (5.15) and (5.26), where tj should 
be taken equal to r n and Shf equal to Sh = b(r n ) — b°(r n ). Calling St = r„ — r° the time delay, we get the following 
relation between Sh and Sh', valid up to 0(St 3 ) corrections : 

Sh = Sh' + b°(r n °) - h°(r n ) = Sh' - Sr^(r n ) + ^Sr^(r° n ) . (5.29) 

Note that St scales typically like \/V in the semi-classical limit. We may thus, to leading order, replace r„ by 
r° and identify Sh with Sh' — St^t(t®) everywhere in the integrand of the right hand side of (5.28), except in 
Y (So(T n ) + Q (i~ n )-Sh) (the combination appearing in the exponential prefactor of (5.15)) which deserves more care. 
It follows from the definition of the action that (again up to 0(St 3 ) corrections) 

S (r B ) = 5 (r°) + St - X -St 2 £°(t°).^(t°) , (5.30) 
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and from (5.29) (together with the condition (r°) ■ Sb 1 = 0) that : 

9\r n ) .Sb = Sr ^(r°) • 5b> - «5t0°(t«) • ^-(r„) + ^ 6° (r° n ) (r° n ) , (5.31) 

When summing these two equations, linear terms in St disappear as expected and we are left after some rearrangements 
with the remainder, of quadratic order in fluctuations : 

1 1 rift® rfVfi rift® 

i {So (r n )+e (r n ).Sb} = + ^ 2 ^(r°) • ^- (r° ) + St ^- (r° ) • Sb , (5.32) 
where we set 5b = Sb' — St — — (r°), so that the time delay St is simply expressed in terms of Sb as St = — <I>2/(t°) Sb. 

(IT 

The end of our theoretical considerations is reached with the following expression for P n (z) : 

P n (^)^e-^°^/ r . = / ^- «S(*;,(t°) -h)dh, (5.33) 

V ^ Vdetf/(rO) 7 (27rr)rf/ 2 V m " ; ; V ; 

where $' 1/ (r) = (b° -b ) 1 / 2 *B*i/(t) (so that *' 1( h = *i/ <*b), W = W + AW, with W defined in the last subsection 
and 

h.AW h = -2($i I • h) (b^b ) 1 / 2 ^- • Bh + • h) 2 ^- • ^- . (5.34) 

We see that the condition of positiveness of the matrix U must be supplemented by the condition of positiveness of 
the restriction W\ of W to the (d — l)-dimensional space orthogonal to &' lh in order to make the instantons found in 
the preceding Section physically meaningful. Provided these two requirements are met, Si(z) is obtained as 

Sl (z) = lim -{Ji(nT) + ilndetZ7(nT) + ^lndetWi(nT)}, (5.35) 
which is the main result of this Section. 

Tracing back all the steps leading to (5.35), one could object to our starting point (5.28) the fact that fluctuations 
of the initial endpoint are not taken into account. This could be done at the expense of rather more cumbersome 
formulae for the reduced path integral Z[h(0) — > h(r)] when both h(0) and h(r) do not vanish. However, we believe 
on physical grounds that the blowing-up associated with the instanton washes out any influence of the fluctuations 
at large scales on the part of the action scaling linearly with the number of steps n. Therefore the expression (5.35) 
should be exact. 



C. Practical implementing and results 

The most difficult part of the computation of s\(z) lies in the evaluation of detf(r) and detW^-r) (as defined 
in the preceding subsection), which requires a good control of all the eigenvalues of these two matrices. Numerical 
instabilities could be avoided for a time long enough to get a precise estimate of -£p\n&etUW\ by using the exact 
expression of the matrix Euler equation (5.22) and the relation (5.25) between W and U. In this problem, there are 
two Goldstone modes associated with uniform rescaling and time translation of the instanton : as a consequence, the 
matrix (b° • b ) 1 / 2 BU(t) (resp. (b° • b ) -1 t B^ 1 WB^ 1 ) is expected to have two eigenvalues scaling like r (resp. like 
1/t) (in order to obtain the simplest transcription of these symmetries, one has to go back to the original fluctuation 
field Sb = (b° • b°) 1 / 2 Bh and the change of variable influences eigenvalues of U and W). When restricting W to 
the (d — l)-dimensional space orthogonal to *' u (= (b° • b ) 1 / 2 ^*!,), one looses one of the eigenvalues scaling like 

1/t, so that there remains an algebraic factor y/r in the product ^AetU{T)\JWi(T), which we had to take away 
by hand in order to make more conspicuous the leading exponential growth of this quantity. To give an idea of the 
accuracy of our procedure and confirm the soundness of the intricate formula (5.35) that was proposed for Si(z), 
we show in Fig. 6 the behaviour of various relevant quantities in the case of model (ii) , and for a moderate scaling 
exponent z = 0.8. It is observed in the picture on the top that the logarithmic derivatives of det U and det W exhibit a 
linear behaviour with almost opposite slopes. It can be shown by considering simpler and exactly solvable models for 
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quadratic fluctuations without inter-shell couplings, that this strange effect mostly reflects the stiff variations suffered 
by the noise variance on any shell, as the latter moves back from the leading edge of the instanton to its rear end. The 
argument is presented in the Appendix C, since it may help the reader to get a feeling for the order of magnitudes at 
play in both U and W matrices. 

It could be tempting at this level to approximate Si (z) as 

lim ^-{X x + I detlnWU}, (5.36) 

r-»oo (IT 2 

This expression would come out if, without great physical justification, it were decided to sum over all positions of 
the final endpoint in the path integral formulation in order to estimate the volume of the basin of attraction of the 
instanton. We found however that the matrix momentum WU develops invariably a negative eigenvalue after some 
time, so that the projection onto the restricted phase space introduced in the previous subsection is a necessary step 
for restoring the statistical stability of the instanton. Further, even before this instability occurs, there was found 
to be a residual r 2 term in IndetWU which forbids any reliable estimate for S\(z) to be deduced from (5.36). The 
picture at the bottom of Fig. 6 shows by contrast that the more precise quantity det U det W\ quickly settles to a 
perfect exponential growth, once algebraic transients have been factored out. Note that the positiveness of both U 
and W\ could be checked in any instance. 

In all the models that we investigated, we found that the first order correction s\ (z) to the action takes an approximate 
parabolic shape of positive or negative concavity, centered around a value of z different from zq ■ In the case of the 
model (ii), the concavity of Si(z) is opposite to the one of Sq(z) and the maximum of Si(z) is reached at a scaling 
exponent Z\ ps 0.6, significantly lower than Zq = 0.72. This means that the minimum of the total action density 
s(z) = sq(z)/T + s\(z) (for values of Y such that s(z) remains concave as it should) is displaced toward the side 
of larger exponents, just as a result of fluctuations. The trend is just the opposite for the model (iii), where s\(z) 
presents this time the same concavity as Sq(z) and a minimum on the left side of Zq. Fig. 7 shows the graph of 
Sq(z) /r + Si(z) that is obtained for the model (ii) and for the value T = 0.58, which we believe to be of some relevance 
for the GOY model (see Section VI). 



VI. COMPARISON WITH NUMERICAL DATA ON THE STATISTICS OF COHERENT EVENTS IN 

THE GOY MODEL 

The systems we have analyzed here were introduced to describe inertial singular structures of the GOY model. To 
check their physical relevance, we first have to define a class of events observed in full simulations of the GOY model 
which are likely to be the best candidates for a description in terms of instantons. It is clear that relative maxima of 
the instantaneous energy flux e n (t) (with e n = fc~ 2 Re{(l — e)Q 2 b n -2b n -ib n + b n -ib n +ib n }) are useful observables for 
tracking the passage of coherent structures across the whole inertial range. But their total number is found to grow 
with n as k n ' , due to the acceleration of time scales typical of the Kolmogorov energy cascade. One may consider that 
they develop on tree- like patterns in the (n,i)-plane, which are renewed at each turn-over of the large scales (totls). 
We say that such trees provide a realization of the propagation of a coherent event from shell no to shell n > no, 
whenever the n — no nodes of the tree closest in time to their supposedly common ancestor on the shell of lower index 
no appear in the order of increasing shell index. In order to discard too weak, and therefore irrelevant, events, we 
imposed that e„ be greater than half the mean energy flux. Figure 8 shows the logarithm of the histogram of the 
logarithmic amplitudes A n = lnle^ 1 / 3 for all relative maxima on one hand and for the restricted class of coherent 
events defined above on the other hand, with no = 5 (far enough from the forcing range) and n = 11 (well in the 
inertial range). The Reynolds number of the simulation is Re = 10 8 . Statistics have been run over 6 x 10 4 totls, and 
in average there are three "coherent lines" for two totls. We note that the statistics of coherent events is very close 
to log-normal for |e„| > 0(1). 

The effective exponent z of a coherent burst after n — no cascade steps is obtained via the relation A n = A no + (n — 
no)(z — 2/3) In 2. If anomalous scaling is preserved in the Re — > oo limit, the pdf of scaling exponents z should behave 
at large cascade lengths as P„(z) ~ e ~ ns ^ z \ where generically s{z) will present a quadratic minimum at some z±, with 
an expansion around z+ that we write as s(z) = a(z — z±) 2 + • • •. The histogram of the variable A n — A no should 
consequently evolve, as n — no increases, towards a Gaussian shape, whose center D n and variance X 2 grow linearly 
with n and relate to z* and a as 

D n ~n In 2(^-2/3), S 2 ~ — (ln2) 2 . (6.1) 

The actual behaviours of D n and — 2E 2 , obtained from a very high Reynolds number simulation {Re = 10 9 , which 
sets the dissipative scale at the shell index = | ln2 Re = 23) are plotted in Fig. 9. Error bars were estimated by 
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varying the range of the quadratic fit to the logarithm of the histogram of A n — A no , as well as the domain of initial 
amplitudes A no used to construct this histogram. It appears that the two graphs are rather far from simple straight 
lines : this is especially true for the variance, whose graph from concave gets convex beyond the shell index 15. From 
the investigation of lower Reynolds numbers Re = 10 8 and Re = 10 7 we could deduce that this transition occurs 
at a shell index n c always of the order of rid — 8 and defines a clear-cut boundary between the inertial range and a 
surprisingly wide pre-viscous range. We believe that the direct action of viscosity on intense bursts starts to show 
up only at the shell index rid — 3, beyond which the local slope of the D„-graph ceases to vary and points to a value 
of the average growth exponent of coherent structures precisely equal to Zq. The fact that n c lies rather far from rid 
means that the cut-off imposed by viscosity exerts a long-range influence on the statistics of the random environment 
seen by a coherent structure. A decent linear regime for both the drift and the variance is observed in the range 
15 < n < 21, from which we get the two estimates a = 29 ± 4 and z± = 0.74 ± 3.10 -3 . Assuming that the fitting range 
10 < n < 16 provides the best clue to the asymptotic scaling of the inertial range, one gets the second set of values 
a = 45 ± 6 and z* = 0.75 ± 3.10~ 3 . It appears that the physics of the previscous range is quite well reproduced within 
our modelling (ii) of the incoherent background. By choosing T = 0.58 (a value hopefully small enough to fall within 
the range of validity of semiclassical approximations), one obtains, as Fig. 7 shows, an almost perfect parabolic shape 
of s(z) = —(Sq(z)/T + Si(z)), with a maximum reached at z* = 0.74 and a curvature a = —2s"(z ic ) = 29. However 
the parameter T cannot be adjusted so as to account for the higher values of a and z* characterizing the inertial 
range. It would seem that in this range of scales it gets necessary to assume some bias in the incoherent fluctuations 
boosting the increase of the renormalized value of z*, while keeping the noise width small. 



VII. CONCLUSION 



We have developed a general scheme for computing numerically self-similar instantons in scale invariant stochastic 
dynamical systems. As concerns the physics of the GOY model, we believe that the bunch of results presented here 
give a strong support to the relevance of an approach focusing from the outset on structures in order to understand 
intermittency and treating the rest of the flow as a noise of weak amplitude. In particular the trend toward log-normal 
statistics of coherent structures is nicely recovered. The detailed study of various stochastic extensions of the GOY 
model shows that the resulting pdf of scaling exponents of singular strucures is very sensitive to the hypothesis made 
on the coupling of noise to the velocity gradient fiels. 

We hope that our approach will be useful for attacking the 3D-Navier-Stokes dynamics along similar lines, once an 
adequate decomposition of the flow into coherent and incoherent parts will have been introduced. An application of 
the method to the Kraignan's model of passive scalar advection formulated on a lattice of shells has already been 
attempted [22]. It has given encouraging results with regard to the validity of a semi-classical analysis, even in 
situations where a small parameter (like T in the present problem) is missing. 
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APPENDIX A 



We carry out in this Appendix the adiabatic approximation alluded to in Section IV B. We look for self-similar 
instantons within the restricted manifold of configurations of the type 

b(r) =e x ^h°(f), (Al) 

where b°(r) is the deterministic solution of scaling exponent zq, f can be thought of as a "proper" time referring 
to the actual position of the pulse. The two variables r(f ) and x(f) parameterize then local changes of speed and 
amplitudes of the pulse, which keeps the same shape as in the absence of noise. Note that if (Al) is to represent a 
self-similar instanton of scaling exponent z ^ Zq, x(f) must obey the constraint 

x(f + T ) - x{f) = (z- z Q ) log Q . (A2) 

We plug now the Ansatz (Al) in the equation of motion (2.11) and project it onto the two directions <&i r = b° and 

^2r = ^jjr. We get by doing so 

df dx T „ ^ ... 

5^ = *"■**' (A3) 
df 

— - 1 = # 2 , ■ fl£ . (A4) 

If the other dimensions of the configuration space are neglected, (A3) and (A4) form a closed two-dimensional system, 
which may be rewritten as 

§ = d, (A5) 
l-g = C 2 , (A6) 

— I/!.- — 

df 

with 



where correlation functions of £j and ( 2 read 



(Ci(r)C j (f l )) = 7 zV ij S(f-f) l<i,j,<2, (A7) 



V ij =*a-B t B* jl . (A8) 
The Gaussian action density associated to one cascade step within this restricted stochastic system is given by 

~ S= 2j df (A9) 
Once expressed in terms of the diffusing variables x(f) and r(f),it becomes 

5 = \ r df{{ % ri + 2 *^ i )i2 + (v-% 2 \ + %(v- i h2 

- [i(V^- 1 )i 2 + (y- 1 ) 22 ]}. (A10) 
The extremization of s with respect to leads to 

^=[(x 2 (V- 1 ) 11 +2x(V- 1 ) 12 + (V- 1 ) 22 )V 22 ] 1/2 , (All) 
and to an effective action for the remaining variable x 

Seff[x(f)}= df{[(x 2 (V-% 1 +2x(V- 1 )y 2 + (V- 1 h 2 )(V- 1 h 2 ] 1/2 

Jo 

- [i(V^- 1 ) 12 + (y- 1 ) 22 ]}. (A12) 

Assuming the coefficients Vn, V\ 2 and V 22 to be almost constant inside the time interval T , one deduces an analytic 
expression for Sq(z) from s e // by just replacing in the integral x by [z — Zq) \ogQ /Tq (which follows from (A2)). One 
gets in particular for large enough z — Zq, 

s (z) ~{z- zo) \ogQ {{{V- 1 )^- 1 )^) 1 ' 2 - {V- 1 )^} , (A13) 

i.e, a linear behaviour as observed for the true solution of model (i). 
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APPENDIX B 



We derive the formal expression of the reduced path integral Z[0 — > h/,r/] given in eq. (5.26) of the Section V of 
the text. The proof is presented in many textbooks on path integrals but, as far as we know, always using a continuous 
definition of time. This leaves some ambiguity in the right equations that matrices U and W should obey, once time 
is discretized for computing purposes. We found that this issue is crucial mostly for evaluating W and preserving its 
symmetry properties. This is why we feel it useful to show how every step of the proof given in the continuum limit 
receives an exact transcription in the discrete time case. 
We start from the quadratic functional : 



S 2 S[h] 



At 

2r 



£ 

i=0 



hj+i - hi 



At 



B 



h,+i + hi 



+ hj+i + hi h i+1 + hi 



(Bl) 



Upon the addition of the boundary term — ^ ^ i=0 {h, + i .Wi + \h i+ i — hj.Wjh,}, it becomes without any approxima- 
tion : 



S 2 S[h] 



At 
2T 



i=N-l 



i=0 



hi 



+1 



At 



Qi\Bi 



Wi + W i+1 . h i+ i + h 



+ 



h i+ i + hj - h i+ i + hj 



-•V, 



(B2) 



where 



Qi = l-^(W i+1 -Wi), 



(B3) 



and 



Vi = V' + t B i B i 



U',. , - Wi 



(% 



^P^)Q-\B, 



Wi + w i+1 , 



(B4) 



We see that <5 2 S[h] reduces to the time integral of a single (positive) square, if and only if Wi is such that for all 
< i < N - 1 



Vi + t B i B i 



W„ 



8+1 



Wi, 



(B5) 



We note that (B5) forces Wi to remain symmetric for all time, provided that Wo (arbitrary at this stage) is chosen 
to be such. In order to solve the above Ricatti matrix equation, one makes the following change of matrix variable 



At 



(B6) 



^From the expression (B3) of Qi, it is easily shown that (B6) is equivalent to the equation (5.25) quoted in the text. 
Furthermore, by multiplying both sides of (B5) by Ui+ ^ i + 1 n the right, one gets for < i < N — 1 



W i+1 U i+1 -WiUi 



u i+ i-u i + {%Bi + v , ) u i + Ui +1 



(B7) 



By half-summing the two relations yielded by (B7) at subsequent values i — 1 and i of the temporal index (with then 
1 < i < N - 1), and using (5.25) after noticing that ^u^-w^u^ = (w^u^+w^-iw^+w^u^) > Qne can 
eliminate W and check that U obeys the matrix Euler equation (5.22), as promised in the text. 



So far, we have proven that, as long as the matrix U may be inverted, the positiveness of S 2 S is guaranteed, since in 
that case the matrix W exists at all times (from (B6)) and allows one to transform the initial quadratic form into the 
time integral of a single square. We show now how these matrices lead to a compact expression of Z[ho = — > h/, r/]. 
We first note that (B6) and (B7) provide 27V relations for 2(7V + 1) unknowns {U , • • • , U N }, {W , • • • , W N }. This 
gives much freedom in the choice of Wq and Uq. In the particular case of h = 0, it is convenient to set Uq = 
and WqUq = 1 (which should be understood as the limit as e — > + of Uq = e and Wo = e _1 , so that Wq is indeed 
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symmetric). A quick inspection of (B6) and (B7) reveals that Ui and Wi behave then respectively as iAr and (2 At) 1 
to leading order in At for 1 < % <C N . One has for instance the exact result W 1 = (At) -1 - B °+' B ° + ('Bo-Bo + V ). 
Recall that the quantity we wish to estimate reads now : 



where we defined the new field ip i+1 (for < i < N — 1) as 

/ l v» a n-i/c , Wi + W i+1 hi + h i+1 
ip i+1 = hi+i - hi - AtQj (B< + ) , 



and the measure of integration Vh as 



N-l 



dh i+ i 



i=0 



(27rrAT) d / 2 



With the help of (B6), the transformation (B9) may be rewritten as 

Ui+i — Ui\ (Ui + Ui+i \ _ 1 



tp i+1 = h i+ i - h, : 
This relation is easily inverted by setting hi = Uid, which gives 

1>i 



(hi + h,+i) . 



(B8) 



(B9) 



(BIO) 



(Bll) 



_f Ui + Ui+i Ui+i -Ui f Ui + Ui+i \ _! Ui+i -Ui \ 
i+1 ~ \ 2 2 V 2 J 2 J 4< 



U i+1 



(d+i Ci) > 



so that we deduce (under the assumption ho = 0) , for < i < N — 1 

hi+i = Ui+i 



3=0 



(B12) 



Since h, + i is linearly related to the V'j+i's °f index j lower than 2, only the diagonal blocks C/j+i I — — 2 i+1 I enter 



the Jacobian of the transformation and one has : 



>N 



dh-i+i 



JV-l 

n 

i=0 



det(Uj + U i+1 ) 
det2tL 



(B13) 



To enforce the boundary condition h/v = h/ at time t/ in terms of the new variables t/ij, we introduce the usual 
integral representation of the (^-function : 



6 d (h f -h N ) 



doc_ ia . [h/ _ E --ii^k^ +i] 
7 (27r) d 



After performing the Gaussian integration over the i/>i's, one arrives at 

JV-l 



Z[0->h /) r / ]= 

8=0 



det([/i + [/i +1 ) ^ 1 



det 2C7i 



^/ det <2, 



e 2r 



7 (2tt)< 



(B14) 



(B1E 



where 



G = U N 



L i=0 



(B16) 
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This awkward non local operator G is greatly simplified when the singular initial condition already mentioned : 
Uq = e, Wq = e -1 with e — > + is adopted. In that case, G is completely dominated by the first term of the series in 
the r.h.s. of (B16) which diverges as e -1 and, to leading order in e, one has 

G ~ ± Un (U^W^U, 1 ) l U N . (B17) 

The summation over a. can then be done and, since G^ 1 vanishes in the e->0 + limit, one gets 

Note that all the manipulations presented in this Appendix were devoid of any approximation. It is finally a 
straightforward matter (details will be skipped here), to check that in the At — > limit, the infinite product in front 

of the exponential in (B18) reduces to ( 2,^ )^ Vdet u ' ma ' sm S thereby (B18) identical to the result (5.26) quoted 
in the text. 



APPENDIX C 

We consider the following quadratic action : 

S2[x n ] = \Y J [ T dr^-, (CI) 
1 „=o ■ / ° a " lrJ 

where, in order to mimic the stochastic models studied in this paper, the variance a n (r) evolves on each shell as 

a„(r) = (b°.b )(^_ 2 ^_ 1 ) 2 , (C2) 

i.e., as the variance of noise, (b° • h°)B nn , for the model (ii). In (C2), b°(r) may be thought of as the deterministic 
self-similar solution, and a n (r) can therefore be cast into the form 

an(r) = e 2A ° T a(r-nT ), (C3) 

where the function a(r) satisfies 

a(r + dT ) = a(r), (C4) 

because of the periodicity of the shell lattice. Let us assume that the instanton is centered around the shell of index 
n = at time r = 0. At its leading edge (n > 0), C n decreases very abruptly as exp —cr n , with r=(>/5 — l)/2 and c 
a constant of order 1. This essential singularity comes from the necessity of balancing with the dominant term 
Q 2 (l — t)C n -iC n -\ of the non-linear kernel of the GOY model in this range of scales. In the trail of the instanton 
(n < 0), one has a much smoother behaviour C n ~ Q nz ° = e nA ° T ° . When the shell lattice is periodized, the leading 
edge and the tail of the instanton have to be glued together and the locus of matching, as well as the residual amplitude 
of C n at that place, will be imposed by the side supporting the slowest variations of C n . We conclude that in a cyclic 
chain containing d shells, most of them reside in the exponential tail of the instanton, so that we may write (again 
under the hypothesis of an instanton initially centered around the origin n = and with a shell index n defined 
between and d — 1) 

a„(0) ~ e i{n - d ^ AoTo , (C5) 

Thus, the range of values spanned by the function a is very large and scales with the total number of shells like e 4dA ° T ° . 
For shells on the exponential ramp, a„(r) first decreases exponentially in time like e~ 2A ° T (because C n decreases like 
e -A r m ^ s re gj on ) anc [ g 0es by a sharp maximum of order e 2nA ° T ° at the time r„ ~ hTq when the center of the 
instanton reaches the corresponding shell. Then it starts again to decrease exponentially. 

Having understood these basic dynamical features, we can compute the matrices U and W introduced in subsection 
VA for the quadratic action given by (CI). To make contact with the normalized field h used in the real problem, 
we switch from the variable x n to the variable y n = x n j ^Ja^. This transforms the original action into 
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d-1 



S2[yn] = \Y, [ dT{y n + \-y n ) 2 . (C6) 
1 „=o Jo 1 an 

Since there is no inter-shell coupling, the matrices U and W are diagonal in the shell index. The extremization of 
S[y„] with respect to y n leads to the couple of first-order differential equations 

Pn=y n + \—y n , (C7) 

Pn = \ —Pn- (C8) 

One has simply [/„„ = y„ and W nn U nn = p n . The solution of Eqs. (C7) and (C8) under the initial conditions 
y n (0) = and p„(0) = 1 is 



Pn( T J - / 77rv ) (C9) 

Va«(0) 

y n (r) = —L=-=L= [ T a n (T')dT'. (CIO) 



\/an(0) \/aJj) Jo 

As long as r < r„ = nTo, the instanton has not passed through the shell of index n and the integral in the right hand 
side of (CIO) is dominated by the neighborhood of the lower bound r = 0. We deduce that for indices n > t/Tq 

Unn(T) = (Cll) 

W nn (r) =2A e- 2A °\ (C12) 

By contrast, when r gets larger than r„ (by some units of time To), the integral in the right hand side of (CIO) is 
dominated by the neighborhood of the time r„ where a n takes its maximal value. We get 



y n {T) ~I\ — j--e 

V a n{0) 



where / is a number of order 1. It follows that for indices n < t/T , 

U nn (r) ~ Ie^ d - n ^ AoTo e A ° T , (C13) 

W nn (r) ~ r 1 e 2nj4oTo e~ 2AoT , (C14) 

where we used the estimate (C5) for a n (0). Since W nn U nn ~ e~ A ° T for n > t/T and ~ e 2dA ° T ° e _ A ° T for n < r/T , we 
conclude that det WU increases exponentially like e dA ° T . But this property is not shared by det U or det W considered 

individually. Indeed, since the number of shells crossed by the instanton increases linearly in time, eqs. (Cll) and 
(C13) show that 



I 



detU ~ ( T7T- ) e 3 ^ dr e-^ ^, (C15) 



while from eqs. (C12) and (C14), 



/ T \ _r l T ° 2 

detW~[—) e -2A dT e +A ^_ (C16) 

This argument captures apparently a good part of the physics of fluctuations around a moving self-similar system, 
though badly treating hybridization effects between neighbouring shells. It also explains how large (resp. small) 
numbers are generated in the spectrum of the matrix U (resp. W) and why in practice one does not have much 
freedom in the choice of the total number of shells d. 
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FIG. 1. On the upper picture (resp. lower) we plot the configurations at a given instant of the right (resp. left) eigenmodes 
in the subspace of maximum Lyapunov exponent Ao around the deterministic solution. 
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FIG. 2. Evolution of the normalized action per unit cascade step Ba(z)/T , as a function of the effective scaling exponent z, 
for the three models studied in this paper. As a guide for the eyes we show the parabola (solid line) "tangent" to the curves 
at the deterministic minimum zo = 0.72. 
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FIG. 3. Configurations of the normalized coherent field C and the random force B£ for three instantons of exponents equal 
to z = 0.75, 0.85, 0.95, obtained with the model (i) (according to the nomenclature defined in the text). Note that the C field 
is only slightly deformed as z increases. 
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FIG. 4. Same as in Fig. 3 but for the model (ii). Note the more pronounced deformation of C upon increasing z, with respect 
to the previous case. 
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FIG. 5. Same as in Fig. 3 but for the model (iii) and only instantons of exponent z = 0.75, 0.85. The physical field now 
changes its sign at the leading edge of the pulse. 
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FIG. 6. Test for the convergence of several relevant quantities entering the calculation of the function Si(z). We show the 
case of the model (ii) at z = 0.8. 
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FIG. 7. Graph (stars) of so{z)/T + si(z) for T = 0.58 and model (ii). The parabolic fit S qua d = a(z — z+) 2 yields z+ = 0.74 
and a = 30. 
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FIG. 8. Histograms of the energy flux in log-log plot, involving all relative maxima of e n or only those associated with 
coherent events. The shell index n = 11, the Reynolds number Re = 10 s and the number of totls is 6 x 10 4 . 
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FIG. 9. The two quantities D n and — 2S^ (encoding the Gaussian central part of the histogram of the growth variable 
A n — A no ) vs n. The Reynolds number Re = 10 9 and the dissipative shell has index n<j = 23. The two pieces of straight line 
show the linear fits that were used to extract the values of and a in the pre-viscous range. 
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